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,  MODELS  FOR  MULTIMEGAWATT  SPACE  POWER  SYSTEMS 


Michael  W .  Edenburn 
Advanced  Power  Systems  Division 
Sandia  National  Laboratories 
Albuquerque,  NM  87185 


ABSTRACT 


This  report  describes  models  for  multimegawatt,  space  power  systems 
which  Sandia' s  Advanced  Power  Systems  Division  has  constructed  to  help 
evaluate  space  power  systems  for  SDI's  Space  Power  Office.  Five  system 
models  and  models  for  associated  components  are  presented  for  both  open 
(power  system  waste  products  are  exhausted  into  space)  and  closed  (no 
waste  products)  systems: 

1.  Open,  burst  mode,  hydrogen  cooled  nuclear  reactor  - 
turboalternator  system; 

2.  Open,  hydrogen- oxygen  combustion  turboalternator  system; 

3.  Closed,  nuclear  reactor  powered  Brayton  cycle  system; 

4.  Closed,  liquid  metal  Rankine  cycle  system;  and 

5.  Closed,  in-core,  reactor  thermionic  system. 

The  models  estimate  performance  and  mass  for  the  components  in  each  of 
these  systems. 
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INTRODUCTION 


Sandia  National  Laboratories'  Advanced  Power  Systems  Division  performs 
system  analysis  for  the  Strategic  Defense  Initiative  (SDI)  Space  Power 
Office's  Independent  Evaluation  Group  (lEG)  to  help  SDI  evaluate 
multimegawatt  space  power  system  concepts.  These  systems  are  ultimately 
planned  to  power  space  based  antiball is tic  missile  weapons  and  their 
associated  surveillance  and  command  satellites .  The  purpose  of  the 
evaluation  is  to  identify  promising  technologies  and  concepts  which  will 
lead  to  more  effective,  less  expensive  space  power  systems.  Identifying 
these  technologies  and  concepts  will  help  the  SDI  direct  their  space 
power  development  program. 

Proposed  multimegawatt  space  power  concepts  take  a  variety  of  forms. 
Primary  energy  sources  proposed  include  reactors,  chemicals,  the  sun, 
and  various  energy  storage  devices  such  as  batteries ,  f l3rwheels ,  and 
large,  cryogenic  inductors.  Turbines  with  generators,  thermionics, 
thermoelectrics,  alkali  metal  thermoelectrics,  thermophotovoltaics , 
magnetohydrodynamics,  and  others  have  been  proposed  as  energy 
converters.  Two  basic  types  of  systems  are  needed- -"continuous"  power 
systems  to  provide  power  over  a  long  time  for  general  service  loads  and 
"burst"  power  systems  to  provide  power  to  weapons  for  a  several  minute 
battle  engagement. 

We  have  constructed  models  for  some  of  the  proposed  power  systems  to 
help  us  in  our  evaluation  efforts.  We  call  these  modeled  systems 
reference  systems  because  we  believe  they  represent  some  of  the  more 
likely  candidates  for  space  power,  and  we  will  use  them  as  a  reference 
of  comparison  to  measure  the  merit  of  other  proposed  system  concepts . 

The  modeled  systems  are  listed  below. 

"Continuous"  power  systems  (Closed) 

Gas  cooled  reactor  powered  Brayton  cycle 
Liquid  metal  cooled  reactor  powered  Rankine  cycle 
Liquid  metal  cooled  reactor  thermionic  system 

"Burst"  power  systems  (Open) 

Gas  cooled  reactor  powered  open  turbine -generator 
Hydrogen -oxygen  combustion  open  turbine -generator 

An  open  system  is  one  that  dumps  turbine  exhaust  into  space.  These 
systems  consist  of  a  power  source,  power  conversion  components,  and 
power  conditioning.  They  do  not  include  the  component  for  which  the 
power  is  being  generated  such  as  a  weapon  or  a  radar.  We  calculate 
system  performance,  system  weight,  and  a  very  rough  estimate  of  system 
cost  using  the  models. 

This  report  will  describe  the  five  models,  and  we  believe  they  will  be 
useful  to  the  SDI's  technology  development  program. 
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GAS  COOLED  REACTOR  POWERED  BURST  SYSTEM 


The  power  source  for  this  system  is  a  hydrogen  cooled  reactor.  Hot 
hydrogen  leaves  the  reactor  and  is  expanded  in  a  turbine  which  produces 
shaft  power  to  run  a  generator  or  alternator.  After  expansion,  the 
hydrogen  is  exhausted  into  space.  Figure  1  shows  a  simple  schematic  of 
this  system.  The  figure  shows  that  the  hydrogen  originates  in  a 
refrigerated  tank  and  that  it  cools  the  weapon,  power  conditioning  unit, 
and  generator  before  entering  the  reactor.  The  weapon,  having  an 
efficiency  less  than  unity,  dumps  heat  into  the  hydrogen  which  enters  as 
a  liquid  and  passes  through  a  supercritical  process  as  it  cools  the 
weapon.  The  amount  of  hydrogen  required  to  cool  the  weapon  is 
determined  by  the  weapon's  cooling  load  and  by  its  prescribed  outlet 
temperature.  If  the  turbine  requires  more  hydrogen  than  the  weapon, 
makeup  hydrogen  is  supplied  from  the  tank  and  bypasses  the  weapon.  If 
the  turbine  requires  less  hydrogen  than  the  weapon,  the  excess  hydrogen 
is  dumped  into  space  after  cooling  the  generator.  The  schematic  also 
shows  a  fl3n*7heel.  This  energy  storage  component  does  not  interact  with 
the  system  in  this  model  but  is  used  to  indicate  that  the  system  will 
store  energy  during  brief  periods  when  the  weapon  is  not  firing  but 
turning  the  power  system  off  is  not  practical.  This  might  happen  when 
the  time  between  bursts  of  shots  is  very  brief  or  when  a  weapon  fault  is 
detected.  The  flyi^heel  contributes  a  small  weight  to  the  system  but  its 
true  size  cannot  be  determined  until  a  dynamic  battle  scenario  model  is 
constructed.  This  type  of  system  is  intended  to  power  an  antiballistic 
missile  weapon.  It  is  believed  that  a  battle  engagement  will  last  less 
than  one-half  hour  and  will  not  be  repeated.  Thus,  the  stored  hydrogen 
must  last  the  duration  of  one  battle  engagement  plus  any  test  time  prior 
to  the  engagement.  The  equations  that  quantify  the  system's  performance 
are  developed  below. 

The  weapon  requires  continuous  electrical  power  (watts) .  If  the 
weapon  pulses,  then  is  the  average  power  needed  but  it  is  still 
supplied  continuously  during  the  battle  engagement.  The  hydrogen  flow 

rate  required  to  cool  the  weapon  is  m^  (kg/sec)  and  is  found  using 
Equation  (1). 


m^[h(T^Q)  -  h(T5tore)]  “  ^w^  ^w^p  > 

where  V'  is  the  weapon  cooling  load  fraction, 

h  is  the  enthalpy  of  hydrogen  (see  Appendix  G) , 

"^wo  is  the  specified  weapon  outlet  temperature  (K) , 

Tgtore  hydrogen  storage  temperature  (20  K,  storage 

pressure  is  assumed  to  be  1  atm) ,  and 
Sp  is  the  pump's  power  per  unit  flow  rate  (Equation  2). 
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EXHAUST 


FIGURE  1.  Hydrogen  Cooled  Reactor, 
Open  Space  Power  System. 
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[Pti(l+dP)  ■  Pstore] 


(2) 


S 


P  " 


0.7  X  Pstore 


where  Pti  is  the  turbine  inlet  pressure  (Pa) . 

dP  is  the  fractional  pressure  drop  across  the  reactor, 

Pstore  is  the  hydrogen  storage  pressure  (Pa), 

Pstore  is  the  hydrogen  storage  density  (kg/m^) ,  and 
0.7  accounts  for  pump  efficiency. 

The  turbine's  outlet  temperature  T^o  depends  on  its  inlet  temperature 
Tti,  its  pressure  ratio  Rp,  and  its  efficiency. 

-a 

Tto  =  Tti  -  Tti(l  -  Rp  )»?t  . 

where  »?t  is  the  turbine's  efficiency,  * 
a  is  R/Cp, 

R  is  hydrogen's  gas  constant,  and 
Cp  is  its  specific  heat  (equation  4) . 

Cp  =  [h(Tti)  ■  k(Tto)]/(Tti  ‘  Tpo)  > 

The  turbine's  flow  rate  is  mj;  (kg/sec),  and  can  be  found  using 
Equation  (5).  * 


where 


mt[h(Tti)  -  h(Tto)] 


Pw  ^  mSp 
ngV-pc  »?g 


(5) 


r?g  is  generator  efficiency,  and 
j/pc  is  power  conditioning  efficiency. 

rii  is  either  the  hydrogen  flow  rate  for  weapon  cooling  or  for 
turbine  power,  whichever  is  greater. 


The  energy  added  to  the  hydrogen  before  it  reaches  the  reactor  is  Qth* 


Qth  “  Pw 


[ 


1 

ngVTpc 


+  ^  -  1 


) 


+ 


The  reactor's  inlet  temperature  is  Tj^£,  which  can  be  found  by 
iteratively  solving  Equation  (7). 


(6) 


*  [  h(T^j) 


-  h(T  ^  ) 

store 


(7) 
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The  reactor's  thermal  power  is  Qj.. 


Q,  =  iiit[h(Tti)  -  h(T,.i)]  .  (8) 

The  parameter  Rp,  the  turbine's  pressure  ratio,  must  be  adjusted  to 
minimize  either  system  weight  or  cost.  As  Rp  increases,  turbine  flow 
rate  is  reduced,  the  reactor's  thermal  power  is  reduced,  and  more 
enthalpy  is  extracted  in  the  turbine.  However,  as  Rp  increases, 
progressively  larger  stages  must  be  added  to  the  turbine  to  accommodate 
the  added  hydrogen  expansion.  Thus,  system  weight,  or  cost,  has  a  "U" 
shaped  curve  when  plotted  against  Rp,  and  there  is  a  minimum  value  which 
must  be  determined. 

The  system  model  consists  of  the  above  equations  and  models  for 
individual  component  weights  and  costs  which  will  be  described  later. 
From  these  equations  and  models,  we  constructed  a  FORTES  program 
(REBRST)  which  runs  on  the  IBM-AT  computer  to  calculate  system 
performance,  weight,  and  cost.  The  model  is  interactive  and  a  listing 
of  the  self  explanatory  input  sequence  is  shown  in  Table  1.  The 
underlined  parameters  are  supplied  by  the  user.  The  model  contains  a 
default  set  of  parameters  that  can  be  replaced  by  following  the  given 
instructions . 

Table  2  shows  a  listing  of  model  output,  The  first  section  of  this 
table  shows  the  optimization  of  pressure  ratio.  As  pressure  ratio 
increases ,  reactor  weight  and  hydrogen  cooling  system  weight  decease  and 
power  conversion  weight  increases  until  pressure  ratio  gets  to  284.9.  A 
pressure  ratio  of  237.4  was  selected  to  minimize  system  weight.  The 
second  section  of  the  table  lists  values  of  selected  parameters  for  the 
optimum  system. 

Table  3  summarizes  the  component  weight  and  cost  models  that  were  used 
in  the  system  model. 

The  reactor  model  was  written  by  Al  Marshall  and  is  described  in 
Marshall  (1986) .  Briefly,  it  calculates  the  fuel  mass  required  for  each 
of  three  limits:  end  of  life  criticality,  burnup,  and  specific  power 
removal.  The  largest  of  the  three  is  selected  and  is  added  to  masses 
for  the  moderator,  structure,  reflector,  pressure  vessel,  shield,  and 
miscellaneous  items  to  get  a  total  reactor  and  shield  mass.  The 
specific  set  of  reactor  parameters  used  in  this  system  model  are  for  a 
hydrogen  cooled,  LiH  moderated,  UC  fuel,  particle  bed,  burst  reactor. 

The  specific  parameter  values  can  be  found  in  Marshall's  report. 

We  have  augmented  Marshall's  model  in  two  areas.  A  specific  power 
calculation  has  been  added,  and  "wraparound"  shield  is  used  when  it  is 
lighter  than  a  planar  shield. 
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Specific  power  is  the  maximum  thermal  power  per  kilogram  of  fuel  that 
can  be  removed  from  the  reactor's  core  without  overheating  the  fuel,  or 
causing  an  excessive  pressure  drop.  The  calculation  of  this  value  is 
described  in  Appendix  A  for  a  hydrogen  or  helium  cooled  particle  bed 
reactor. 

The  wraparound  shield  model  is  explained  in  Appendix  B.  Marshall's 
model  calculates  the  weight  of  a  disc  shaped  shadow  shield.  When  the 
protection  cone  angle  is  large,  a  wraparound  shield  will  be  lighter  than 
a  disc  shadow  shield  and  is  substituted  when  this  is  the  case. 

The  turbine  model  is  described  in  Appendix  C.  It  is  based  on  a  stage- 
by-stage  size,  weight,  and  efficiency  computation  with  temperature 
dependent  material  strengths,  working  fluid  properties,  and  working 
fluid  energy  extraction.  Table  3  describes  the  models  for  the  fl^heel 
and  for  the  generator  which  is  assumed  to  have  a  95  percent  efficiency. 
Power  conditioning  weight  is  assumed  to  be  0.2  kg/kW.  We  estimate  that 
beam  weapons  will  require  power  conditioners  that  weigh  close  to 
0  5  kg/kW  and  that  kinetic  energy  weapons  will  require  almost  no  power 
conditioning  weight.  Thus,  0.2  kg/kW  is  used  to  represent  the  range  of 
estimates. 

The  cooling  system  consists  of  hydrogen,  an  insulated  hydrogen  tank,  a 
refrigeration  system  to  keep  the  hydrogen  at  cryogenic  temperatures,  and 
a  meteoroid  shield.  The  model  for  this  system  is  described  in 
Appendix  D.  The  turbine  and  the  weapon  may  require  different  quantities 
of  hydrogen.  The  model  uses  the  greater  of  the  two  values  to  calculate 
cooling  system  weight. 

The  weapon  is  simply  modeled  as  a  cooling  load.  If  its  cooling  load  is 
50  percent,  then  half  of  the  energy  used  by  the  weapon  must  be  extracted 

as  heat. 

The  model  also  includes  cost.  While  our  estimates  for  these  are  not  yet 
very  accurate,  they  have  been  included  because  we  recognize  that  cost 
will  become  the  most  important  factor  by  which  to  judge  a  system.  The 
major  cost  for  previously  launched  small  power  system  was  launch  cost 
which  is  proportional  to  weight,  but  for  multimegawatt  size,  lower 
specific  weight  systems,  the  production  cost  of  the  system  is  also 
important.  Our  estimates  for  launch  cost  are  based  on  present  shuttle 
launch  cost.  Aviation  Week  (1985)  estimates  that  a  commercially 
competitive  shuttle  launch  would  cost  $128  million.  If  a  heavy  lift 
shuttle  could  lift  100  metric  tons  into  low  earth  orbit,  then  the  launch 
cost  would  be  $1280/kg.  This  is  the  launch  cost  used  by  the  model. 
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TABLE  1 


THIS  PROGRAM  MODELS  AN  OPEN  CYCLE  SPACE  POWER 
SYSTEM,  THE  SYSTEM  CONSISTS  OF  A  HYDROGEN  COOLED 
REACTOR,  A  TURBINE,  A  FLYWHEEL,  A  GENERATOR, 

A  POWER  CONDITIONER,  A  WEAPON  SUBSYSTEM, 

AND  A  HYDROGEN  STORAGE  SUBSYSTEM. 

PRESSURE  RATIO  IS  OPTIMIZED  TO  GET  EITHER 
MINIMUM  SYSTEM  WEIGHT  OR  COST. 

THE  WEAPON  IS  COOLED  BY  SUPERCRITICAL  HYDROGEN 
STORED  AT  1  ATM.  AS  A  LIQUID  AND  PUMPED  TO  REACH 
TURBINE  PRESSURE. 

THE  COOLANT  ALSO  COOLS  THE  POWER  CONDITIONER 
AND  THE  GENERATOR  BEFORE  ENTERING  THE  REACTOR. 


ENTER  VALUES  FOR  WEAPON  INPUT  ELECTRICAL  POWER 
IN  MW  AND  OPERATING  TIME  IN  HOURS.  100  .2 

ENTER  "WEIGHT”  OR  "COST"  AS  PARAMETER  TO  BE 
MINIMIZED  (IN  CAPITAL  LETTERS) .  WEIGHT 

THE  FOLLOWING  ARE  DEFAULT  PARAMETER  VALUES, 
THEY  MAY  BE  CHANGED  IF  YOU  WISH, 


CYCLE  PARAMETERS: 

1.  1200.00000000 

2.  4.00000000 

3.  l.OOOOOOOOE+06 

4.  0.94999999 

5.  0.94999999 

6.  0.50000000 

7.  300.00000000 

8.  10.00000000 


TURBINE  INLET  TEMP,  K 
TURBINE  INLET  PRESSURE,  MPa 
MAXIMUM  PRESSURE  RATIO 
GENERATOR  EFFICIENCY 
POWER  CONDITIONING  EFFICIENCY 
WEAPON  COOLING  LOAD  FRACTION 
WEAPON  OUTLET  TEMP,  K 
FLYWHEEL  ENERGY  STORAGE  TIME,  SEC 


TURBINE  PARAMETERS: 


9. 


4.00000000  NUMBER  OF  TURBINES 


10, 

1.00000000 

2-CARB 

11. 

1200.00000000 

12. 

1200.00000000 

13. 

10000.00000000 

14. 

5.00000000 

REACTOR  PARAMETERS: 

15.  0.93000001 

16.  1.00000000 

17.  0.40000001 

18.  2.19999999E-02 

19.  0.10000000 

20.  0.25000000 

21.  65.50000000 

22.  93.00000000 

23.  5610.00000000 


FRACTIONAL  FUEL  ENRICHMENT 
CRITICAL  COMPACT  MASS,  Kg 
FUEL  +  MODERATOR  VOL  FRACT 
FUEL  BED  LENGTH,  m 
REACTOR  PRSUR  DROP,  FRACT  OF  PIN 
FUEL  BURNUP  FRACTION  LIMIT 
MODERATOR-TO-FUEL  RATIO 
MODERATOR  MOLECULAR  WEIGHT 
MODERATOR  DENSITY,  Kg/itl3 
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TABLE  1  (cont.) 


24.  4.99999992E+16 

25.  25.00000000 

26.  15.00000000 

27.  2.00000000 

28.  1. OOOOOOOOE+07 


ALLOWED  PAYLOAD  NEUTRON  DOSE,  nvt  . 
PAYLOAD  SEPARATION  DISTANCE,  m 
PROTECTION  CONE  HALF  ANGLE,  DEG 
NEUTRON  SHIELD  MATL-B4C=1,  LIH=2 
ALLOWED  PAYLOAD  GAMMA  DOSE,  R 


ENTER  THE  NUMBER  OF  PARAMETERS  YOU  WISH  TO  CHANGE.  _3_ 
ENTER  THE  3  PARAMETER  NUMBERS.  1  11  12 

ENTER  THE  3  PARAMETER  VALUES.  1150  900  1150 


CYCLE  PARAMETERS: 

1.  1150.00000000 

2.  4.00000000 

3.  l.OOOOOOOOE+06 

4.  0.94999999 

5.  0.94999999 

6.  0.50000000 

7.  300.00000000 

8.  10.00000000 


TURBINE  INLET  TEMP,  K 
TURBINE  INLET  PRESSURE,  MPa 
MAXIMUM  PRESSURE  RATIO 
GENERATOR  EFFICIENCY 
POWER  CONDITIONING  EFFICIENCY 
WEAPON  COOLING  LOAD  FRACTION 
WEAPON  OUTLET  TEMP,  K 
FLYWHEEL  ENERGY  STORAGE  TIME,  SEC 


TURBINE  PARAMETERS: 


9. 


4.00000000  NUMBER  OF  TURBINES 


10, 

1.00000000 

2-CARB( 

11. 

900.00000000 

12. 

1150.00000000 

13. 

10000.00000000 

14. 

5.00000000 

1-SUPERALLOY , 


REACTOR  PARAMETERS: 

15,  0.93000001 

16.  1.00000000 

17.  0.40000001 

18.  2.19999999E-02 

19.  0.10000000 

20.  0.25000000 

21.  65.50000000 

22.  93.00000000 

23.  5610.00000000 

24.  4.99999992E+16 

25.  25.00000000 

26.  15.00000000 

27.  2.00000000 

28.  l.OOOOOOOOE+07 


FRACTIONAL  FUEL  ENRICHMENT 
CRITICAL  COMPACT  MASS,  Kg 
FUEL  +  MODERATOR  VOL  FRACT 
FUEL  BED  LENGTH,  m 
REACTOR  PRSUR  DROP,  FRACT  OF  PIN 
FUEL  BURNUP  FRACTION  LIMIT 
MODERATOR-TO-FUEL  RATIO 
MODERATOR  MOLECULAR  WEIGHT 
MODERATOR  DENSITY,  Kg/m3 
ALLOWED  PAYLOAD  NEUTRON  DOSE,  nvt 
PAYLOAD  SEPARATION  DISTANCE,  m 
PROTECTION  CONE  HALF  ANGLE,  DEG 
NEUTRON  SHIELD  MATL-B4C=1,  LIH=2 
ALLOWED  PAYLOAD  GAMMA  DOSE,  R 


ENTER  THE  NUMBER  OF  PARAMETERS  YOU  WISH  TO  CHANGE. 
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TABLE  2 


PRESSURE 

RCT+SHLD 

RATIO 

WT  (Kg) 

1.200 

74782. 

1.440 

36241. 

1.728 

24741. 

2.074 

19177. 

2.488 

15880. 

2.986 

13693. 

3.583 

12132. 

4.300 

10959. 

5.160 

10044. 

6.192 

9309. 

7.430 

8705. 

8.916 

8200. 

10.699 

7770. 

12.839 

7399. 

15.407 

7076. 

18.488 

6792. 

22.186 

6541. 

26.623 

6316. 

31.948 

6113. 

38.338 

5931. 

46.005 

5764. 

55.206 

5613. 

66.247 

5496. 

79.497 

5391. 

95.396 

5294. 

114.476 

5205. 

137.371 

5122. 

164.845 

5046. 

197.814 

4975. 

237.377 

4902. 

284.852 

4814. 

237.377 

4902. 

POW 

CONV 

COOL  SYS 

WT 

(Kg) 

WT  (Kg) 

38310. 

284803. 

37121. 

136823. 

36836. 

92074. 

36750. 

70299. 

36741. 

57375. 

36773. 

48800. 

36828. 

42688. 

36900. 

38110. 

36984. 

34550. 

37078. 

31703. 

37180. 

29374. 

37290. 

27434. 

37407. 

25793. 

37530. 

24388. 

37660. 

23171. 

37797. 

22108. 

37939. 

21170. 

38089. 

20339. 

38244. 

19596. 

38406. 

18928. 

38575. 

18326. 

38751. 

17780. 

38934. 

17283. 

39124. 

16828. 

39321. 

16411. 

39527. 

16027. 

39740. 

15673. 

39961. 

15346. 

40191. 

15042. 

40431. 

14831. 

40781. 

14831. 

40431. 

14831. 

TOTAL 

TOTAL 

WT  (Kg) 

COST  (M$) 

437685. 

1142. 

231203. 

701. 

169016. 

569. 

138848. 

505. 

120996. 

467. 

109192. 

442. 

100813. 

425. 

94565. 

412. 

89736. 

401. 

85899. 

393. 

82785. 

387. 

80216. 

382. 

78067. 

377. 

76249. 

373. 

74698. 

370. 

73366. 

367. 

72215. 

365. 

71217. 

363. 

70349. 

361. 

69592. 

360. 

68932. 

359. 

68358. 

358. 

67884. 

357. 

67477. 

356. 

67129. 

356. 

66835. 

355. 

66589. 

355. 

66389.- 

355. 

66229. 

355. 

66179. 

355. 

66469. 

356. 

66179. 

355. 

10.73518658 

10.68002987 

1150.00000000 

451.15634155 

0.92879885 

369.53829956 

237.37660217 

124.45237732 


1.16099336E-03 

24.81669044 

24.81677246 

24.81793213 

3.37396975E+06 

47.95184708 


THERMODYNAMIC  PARAMETERS 
WEAPON  FLOW  RATE-Kg/S 
TURBINE  FLOW  RATE-Kg/s 
TURBINE  INLET  TEMP-K 
TURBINE  OUTLET  TEMP-K 
PUMP  POWER-MW 
REACTOR  INLET  TEMP-K 
PRESSURE  RATIO 
REACTOR  THERMAL  POWER-MW 


REACTOR  PARAMETERS 
BURNUP  MASS-Kg. 

INITIAL  CRITICAL  MASS-Kg 

END  CRITICAL  MASS-Kg 

TOTAL  BRNUP+CRIT  MASS-Kg 

SPECIFIC  POWER-W/Kg 

MASS  FOR  SPECIFIC  POWER  LIM-Kg 
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TABLE  2  (cont.) 


MASS  FOR  ALLOWED  BURNUP-Kg 
FUEL  MASS-Kg 
MODERATOR  MASS-Kg 
STRUCTURE  MASS-Kg 
REFLECTOR  MASS-Kg 
PRESSURE  VESSEL  MASS-Kg 
MISCELLANEOUS  MASS-Kg 
TOTAL  REACTOR  MASS-Kg 
NEUTRON  SHIELD  THICKNESS-m 
NEUTRON  SHIELD  MASS-Kg 
GAMMA  SHIELD  THICKNESS-m 
GAMMA  SHIELD  MASS-Kg 
TOTAL  SHIELD  MASS-Kg 


5.61456382E-03 

47.95184708 

1155.96496582 

1373.03930664 

871.28448486 

249.86604309 

1203.91687012 

4902.02343750 

O.OOOOOOOOE-01 

O.OOOOOOOOE-01 

O.OOOOOOOOE-01 

O.OOOOOOOOE-01 

O.OOOOOOOOE-01 


HYDROGEN  STORAGE  PARAMETERS 
TURBINE 
TANK  VOLDME-m3 

REFRIGERATION  SYSTEM  POWER-MW 
HYDROGEN  WEIGHT-Kg 
TANK  WEIGHT-Kg 
INSULATION  WEIGHT-kg 
REFRIG  SYST  WEIGHT-KG 
METEOROID  SHIELD  WEIGHT-kg 
WEAPON 

TANK  VOLUME-m3 

REFRIGERATION  SYSTEM  POWER-MW 
HYDROGEN  WEIGHT-Kg 
TANK  WEIGHT-Kg 
INSULATION  WEIGHT-kg 
REFRIG  SYST  WEIGHT-KG 
METEOROID  SHIELD  WEIGHT-kg 


108.30452728 

7.05457712E-03 

7689.62158203 

219.85820007 

352.72885132 

484.09591675 

6014.22802734 

108.86386871 

7.07885763E-03 

7729.33447266 

220.99363708 

353.94287109 

485.36471558 

6040.93066406 


POWER  CONVERSION  PARAMETERS 
TURBINE  STAGES 
TURBINE  SPEED-RPM 
TURBINE  EFFICIENCY 
TURBINE  PRESSURE  RATIO 
TURBINE  WEIGHT-Kg 
GENERATOR  WEIGHT-Kg 
FLYWHEEL  WEIGHT-Kg 
POWER  CONDITIONING  WEIGHT-Kg 


19 

10000.00000000 

0.78069597 

237.37660217 

6733.51806641 

10619.19628906 

3077.87011719 

20000.00000000 
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TABLE  2  (cont.) 


WEIGHT  AND  COST  SUMMARY 
REACTOR  WEIGHT-Kg 
SHIELD  WEIGHT-Kg 
POWER  CONVERSION  WEIGHT-Kg 
COOLING  SUBSYSTEM  WEIGHT-Kg 
MISC.  WEIGHT-Kg 
TOTAL  WEIGHT-Kg 


4902-02343750 

O.OOOOOOOOE-01 

40430.58593750 

14830.56640625 

6016.31787109 

66179.49218750 


REACTOR+SHIELD  COST-M$ 
POWER  CONVERSION  COST-M$ 
COOLING  SUBSYSTEM  COST-M$ 
FLYWHEEL  COST-M$ 


21.56890297 

245.32461548 

0.31664482 

3.07787037 


LAUNCH  COST-M$  ,=  84.70974731 

TOTAL  COST-M$  =  354.99777222 

Execution  terminated  :  0 
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TABLE  3 

Reactor  Powered  Burst  System 
Weight  and  Cost  Algorithm  Simmnary 


Component 

Reactor  and  Shield 

weient 

Marshall ' s  algorithm 
(Marshall,  1986) 

$4400/kg 

(estimate) 

Turbine 

See  Appendix  C 

$ 2000/kg 
(Gerry,  1985) 

Generator 

0.1  kgAW  (Gerry,  1985) 

$ 3000/kg 
(Gerry,  1985) 

Power  Conditioning 

0.2  kg/kW  (estimate) 

$10, 000/kg 
(estimate) 

Hydrogen  Cooling 
Subsystem 

Wh  =  mT  i 

V  =  Wh/^ 

Wt  =  2.03  V 

Wi  =  15.5  VO-667 

Wj.  =  12.6  VO -667+109x10** 

[0.0468(log(2.78V**0.667))**2.9] 

W3  =  107  VO -86 

Wh  =  hydrogen  mass,  kg 

m  =  hydrogen  flow  rate,  kg/s 

T  =  operating  time ,  s 

V  =  volume  of  hydrogen,  m^ 

p  =  density  of  liquid  H2,  71  kg/m® 

Wt  =  tank  weight,  kg 

Wi  =  insulation  weight,  kg 

Wr  “  refrigeration  weight,  kg 

Wg  =  meteoroid  shield  weight,  kg 

See  Appendix  D 

... 

hydrogen  $2 8/kg 
(Bents,  1984) 
tank  $20/kg 
(estimate) 
refrigeration 
$640/kg  (est.) 

insulation 
$ 100/kg  (est.) 

shield  $10/kg 
(estimate) 

Miscellaneous 

10%  of  subtotal 

Launch  cost 

! 

$1280/kg 
(Aviation  Week, 
1985,  and  heavy 
lift  shuttle) 
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H2-02  COMBUSTION  POWERED  BURST  SYSTEM 


This  system  is  similar  to  the  reactor  powered  burst  system,  but  instead 
of  a  reactor,  it  has  an  oxygen  subsystem  and  a  hydrogen- oxygen  combustor 
(see  Figure  2) .  The  turbine  in  this  system  uses  the  combustion 
products,  a  mixture  of  steam  and  hydrogen,  as  a  working  fluid.  If  a 
stoichiometric  mixture  of  hydrogen  and  oxygen  is  burned,  the  combustion 
product  temperature  will  be  too  high  for  a  turbine  inlet  temperature. 

To  cool  the  combustion  products,  excess  hydrogen  is  used  in  the 
combustion  mixture.  The  ratio,  R,  of  hydrogen  to  oxygen  determines  the 
combustion  product  temperature  (which  is  also  the  turbine  inlet 
temperature  .  Tti).  Thus,  for  a  desired  combustion  product  temperature, 
there  is  a  required  ratio  of  hydrogen  to  oxygen  as  specified  by 
equation  9  and  derived  in  Appendix  E. 


0 . 125  [hH(T)  -had,,)  ]+l .  125hcoMB-ho(T^)+ho(T3o)  -1.125  [hs(T)  -hs(T J  J+Sp^ 

R - 

had)  -  hadia) 


(9) 


where  hcoMB  enthalpy  of  combustion  (13400  J/kg  H2O)  , 

Tih  is  the  temperature  of  the  hydrogen  entering  the  combustion 
chamber, 

T  is  the  combustion  product  temperature, 
ho  is  the  enthalpy  of  oxygen, 
hs  is  the  enthalpy  of  steam, 
ha  is  the  enthalpy  of  hydrogen, 

Spo  is  the  oxygen  pump  power  per  unit  of  oxygen  flow, 

Tgo  is  the  oxygen's  storage  temperature,  and 
Tj.  is  the  reaction  temperature  (300  K) . 


Enthalpy  values  are  specified  in  Appendix  G. 

The  thermodynamic  properties  for  the  combustion  products  are  specified 
in  Appendix  E.  They  are  used  to  calculate  turbine  performances  and 
weight. 


Hydrogen  is  used  to  cool  the  weapon,  power  conditioning  unit,  and 
generator  before  it  enters  the  combustor.  The  amount  of  hydrogen  needed 
depends  on  weapon  efficiency  and  weapon  outlet  temperature.  The  amount 
of  hydrogen  needed  by  the  turbine  depends  on  its  inlet  temperature, 
pressure  ratio,  efficiency,  and  on  the  ratio  of  hydrogen  to  steam  in  the 
working  fluid.  It  may  need  either  more  or  less  hydrogen  than  the 
weapon.  If  it  needs  more,  extra  hydrogen  will  be  supplied  by  the 
hydrogen  storage  subsystem.  If  it  needs  less,  some  weapon  coolant 
hydrogen  will  be  dumped  after  cooling  the  weapon,  power  conditioning 
unit,  and  generator.  The  hydrogen  needed  on  the  platform  is  the  greater 
of  that  needed  by  the  weapon  and  that  needed  by  the  turbine. 
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EXHAUST 


FIGURE  2.  H2-O2  Combustion, 

Open  Space  Power  System. 
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The  system  includes  a  flywheel  to  accommodate  short  energy  transients. 

It  has  no  effect  on  system  performance  in  this  model.  It  does,  however, 
add  a  small  amount  of  weight  to  the  system. 

The  weapon  is  presumed  to  require  continuous  electrical  power,  P^,  to 
fight  a  battle  over  a  several -minute  time  interval.  The  flow  rate  of 
coolant  hydrogen  needed  by  the  weapon  is  specified  by  Equation  (10) . 


^hwI^H^'^wo)  "  ^hC^sh)]  “  ^w^  ^hw^ph  >  (10) 

where  ih^w  Is  the  mass  flow  rate  of  hydrogen  to  the  weapon  (kg/sec) , 

^  is  the  weapon's  cooling  load  fraction, 
turbine  power,  whichever  is  greater, 

^wo  is  the  specified  weapon  outlet  temperature  (K) , 

Tgh  Is  the  temperature  at  which  liquid  hydrogen  is  stored 
(20  K) ,  and 

Sph  is  the  power  divided  by  flow  rate  for  the  hydrogen  pump . 

The  turbine's  outlet  temperature  T^o  depends  on  its  inlet  temperature 
T^i,  its  pressure  ratio  Rp,  and  its  efficiency 

Ito  ^  Iti  “  1ti(l  ‘  ^p  )^t  >  (11) 

where  oc  is  Rg/Cp  where  Rg  is  the  gas  constant  and  Cp  is  the  specific 
heat  for  the  steam-hydrogen  mixture.  Cp  is  calculated  using 
weighted  enthalpy  change  values  (see  Appendix  G) . 

The  turbine's  total  flow  rate  is  m^. 

*t[hc(Tcl)  -  hc(Tto)l  -  ^  ^  ,  (12) 

^g^pc  ^g 

where  h^  is  the  enthalpy  of  the  hydrogen- steam  mixture  (Appendix  E) , 

riiQ  is  the  mass  flow  rate  of  oxygen, 

r^g  is  the  generator  efficiency, 

r/pc  is  the  power  conditioning  efficiency,  and 

mh  is  the  flow  rate  of  hydrogen  for  either  weapon  cooling  or 
turbine  power,  whichever  is  greater. 

(^ti“^s) 

Sp^  =  j  (13) 

0.7pH 


(Pti'Ps) 

Spo  “  ,  (1^) 

0.7Po 
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where 


Pti  is  the  turbine  inlet  pressure  (Pa), 

Pg  is  the  storage  pressure  (100,000  Pa), 

PH  is  the  density  of  liquid  hydrogen  (71  kg/m^) , 

PO  is  the  density  of  liquid  oxygen  (1142  kg/m^) ,  and 
0.7  represents  pump  efficiency. 


The  turbine's  hydrogen  and  oxygen  flow  rates  are  found  using  Equations 
(15)  and  (16),  respectively. 


mht  -  lilt  i+R 


(15) 


nin  "" 


mt 


1+R 


(16) 


The  enthalpy  change  of  the  coolant  hydrogen  before  it  enters  the 
combustor  is  given  by  Equation  (17) . 


Qth 


Pw  f  — ^  +  V  -  1  1  +  +  ri-oSpo  f  ^  -  1  1  •  (17) 

L  ^gVipc  )  I  'jg  J 


The  hydrogen  enters  the  combustor  at  temperature  Tj-^  which  can  be  found 
by  iterating  Equation  (18) . 


*h  [  liH(Tri)  -  hH(Tsh)  =  Qth  ] 


(18) 


Recall  the  pressure  ratio  parameter,  Rp,  that  was  used  to  determine 
turbine  outlet  temperature.  This  parameter  is  a  variable  and  must  be 
adjusted  to  minimize  either  system  weight  or  cost.  When  Rp  increases, 
turbine  enthalpy  extraction  increases,  and  the  flow  rates  of  hydrogen 
and  oxygen  decrease.  However,  when  Rp  increases,  the  turbine  becomes 
heavier  because  it  must  have  additional  stages.  Because  of  this,  there 
is  an  optimum  value  of  Rp  that  minimizes  either  system  weight  or  system 
cost. 

The  above  equations,  along  with  weight  and  cost  models  (discussed 
later) ,  for  each  of  the  system  components  have  been  combined  into  a 
system  computer  model  that  calculates  performance,  weight,  and  cost. 

The  system  model  (HOBRST)  is  written  in  FORTRAN  and  runs  on  an  IBM-AT 
computer.  It  is  interactive  and  a  listing  of  the  self  explanatory  input 
sequence  is  shown  in  Table  4.  The  underlined  parameters  are  supplied  by 
the  user.  The  model  contains  a  default  set  of  parameters  that  can  be 
replaced  by  following  the  instructions  given  while  running  the  program. 
Table  5  is  a  listing  of  model  output.  The  first  part  of  the  table  shows 
how  the  pressure  ratio  is  optimized.  The  optimum  value  selected  was 
114.  The  rest  of  the  table  summarizes  values  of  selected  parameters 
from  the  optimum  system. 
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The  weight  and  cost  models  for  this  system  are  sximmarized  in  Table  6  and 
consist  of  hydrogen  and  oxygen  subsystems,  a  turbine,  a  generator,  power 
conditioning,  and  a  flywheel.  The  turbine  model  is  described  in  detail 
in  Appendix  C.  The  generator  and  power  conditioning  unit  are  both  95 
percent  efficient.  The  cryogen  (hydrogen  and  oxygen)  storage  subsystems 
consist  of  the  stored  liquid,  tanks,  multilayer  insulation, 
refrigeration  units,  and  meteoroid  shields.  The  weights  of  hydrogen  and 
oxygen  are  calculated  by  multiplying  their  flow  rates  by  the  system's 
operating  time.  The  other  cryogen  storage  component  weights  are 
addressed  in  detail  in  Appendix  D. 

The  weapon  is  modeled  as  shown  in  equation  10.  It  contributes  no 
weight,  other  than  its  required  coolant,  to  the  system. 

The  model  also  includes  cost.  The  cost  parameters  in  the  model  are  very 
crude  at  present,  but  they  represent  an  attempt  to  consider  the 
parameter  that  should  be  the  most  important  discriminator  in  the  future. 
The  launch  cost  used  in  the  model  is  a  reasonable  estimate,  $1280/kg, 
based  on  current  shuttle  experience  (Ref.  5)  projected  to  a  heavy  lift 
vehicle  that  can  put  100  metric  tons  in  orbit. 
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TABLE  4 


THIS  PROGRAM  MODELS  AN  OPEN  CYCLE  SPACE  POWER 
SYSTEM.  THE  SYSTEM  CONSISTS  OF  A  HYDROGEN — OXYGEN 
COMBUSTION  TURBINE,  A  FLYWHEEL,  A  GENERATOR, 

A  POWER  CONDITIONER,  A  WEAPON  SUBSYSTEM, 

AND  HYDROGEN  AND  OXYGEN  STORAGE  SUBSYSTEMS. 
PRESSURE  RATIO  IS  OPTIMIZED  TO  GET  EITHER 
MINIMUM  SYSTEM  WEIGHT  OR  COST. 

THE  WEAPON  IS  COOLED  BY  SUPERCRITICAL  HYDROGEN 
STORED  AT  1  ATM.  AS  A  LIQUID  AND  PUMPED  TO  REACH 
TURBINE  PRESSURE. 

THE  HYDROGEN  COOLS  THE  POWER  CONDITIONER 

AND  THE  GENERATOR  BEFORE  ENTERING  THE  COMBUSTOR. 


ENTER  VALUES  FOR  WEAPON  INPUT  ELECTRICAL  POWER 
IN  MW  AND  OPERATING  TIME  IN  HOURS.  100  >2 

ENTER  "WEIGHT”  OR  "COST”  AS  PARAMETER  TO  BE 
MINIMIZED  (IN  CAPITAL  LETTERS) .  WEIGHT 

THE  FOLLOWING  ARE  DEFAULT  PARAMETER  VALUES. 
THEY  MAY  BE  CHANGED  IF  YOU  WISH. 


CYCLE  PARAMETERS; 

1.  1200.00000000 

2.  4.00000000 

3.  l.OOOOOOOOE+06 

4.  0.94999999 

5.  0.94999999 

6.  0.50000000 

7.  300.00000000 

8.  10.00000000 


TURBINE  INLET  TEMP,  K 
TURBINE  INLET  PRESSURE,  MPa 
MAXIMIUM  PRESSURE  RATIO 
GENERATOR  EFFICIENCY 
POWER  CONDITIONING  EFFICIENCY 
WEAPON  COOLING  LOAD  FRACTION 
WEAPON  OUTLET  TEMP,  K 
FLYWHEEL  ENERGY  STORAGE  TIME,  SEC 


TURBINE  PARAMETERS: 

9.  4.00000000 

10.  1.00000000 


NUMBER  OF  TURBINES 
TURBINE  MATERIAL;  1 -SUPERALLOY, 
2 -CARBON  COMPOSITE 


11. 

12. 

13. 

14. 


1200.00000000 

1200.00000000 

10000.00000000 

4.00000000 


MAXIMUM  DISK  TEMPERATURE,  K 
MAXIMUM  BLADE  TEMPERATURE,  K 
MAXIMUM  TURBINE  SPEED,  RPM 
WORK  COEFFICIENT 


ENTER 

THE  NUMBER  OF 

PARAMETERS  YOU  WISH 

TO  CHANGE. 

3 

ENTER 

THE 

3 

PARAMETER 

NUMBERS. 

1  11  12 

ENTER 

THE 

3 

PARAMETER 

VALUES . 

1150  900 

1150 
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TABLE  4  (cont.) 


CYCLE  PARAMETERS: 

1.  1150.00000000  TURBINE  INLET  TEMP,  K 

2.  4.00000000  TURBINE  INLET  PRESSURE,  MPa 

3.  l.OOOOOOOOE+06  MAXIMIUM  PRESSURE  RATIO 

4.  0.94999999  GENERATOR  EFFICIENCY 

5.  0.94999999  POWER  CONDITIONING  EFFICIENCY 

6.  0.50000000  WEAPON  COOLING  LOAD  FRACTION 

7.  300.00000000  WEAPON  OUTLET  TEMP,  K 

8.  10.00000000  FLYWHEEL  ENERGY  STORAGE  TIME,  SEC 


TURBINE  PARAMETERS: 

9.  4.00000000  NUMBER  OF  TURBINES 

10.  1.00000000  TURBINE  MATERIAL;  l-SUPERALLOY, 

2 -CARBON  COMPOSITE 

11.  900.00000000  MAXIMUM  DISK  TEMPERATURE,  K 

12.  1150.00000000  MAXIMUM  BLADE  TEMPERATURE,  K 

13.  10000.00000000  MAXIMUM  TURBINE  SPEED,  RPM 

14.  4.00000000  WORK  COEFFICIENT 


ENTER  THE  NUMBER  OF  PARAMETERS  YOU  WISH  TO  CHANGE.  0 


-19- 


TABLE  5 


PRESSURE 

02-SYSTEM 

H2 -SYSTEM 

RATIO 

WT  (Kg) 

WT  (Kg) 

1.200 

194770. 

244903. 

1.440 

90008. 

119658. 

1.728 

58537. 

81011. 

2.074 

43369. 

62078. 

2.488 

34443. 

50804. 

2.986 

28565. 

43310. 

3.583 

24404. 

37964. 

4.300 

21304. 

33956. 

5.160 

18908. 

30840. 

6.192 

17000. 

28348. 

7.430 

15447. 

26310. 

8.916 

14159. 

24614. 

10.699 

13074. 

23180. 

12.839 

12148. 

21952. 

15.407 

11349. 

20890. 

18.488 

10653. 

19963. 

22.186 

10042. 

19147. 

26.623 

9502. 

18423. 

31.948 

9014. 

17778. 

38.338 

8590. 

17198. 

46.005 

8202. 

16676. 

55.206 

7851. 

16204. 

66.247 

7532. 

15774. 

79.497 

7242. 

15381. 

95.396 

6976. 

15022. 

114.476 

6764. 

14807. 

137.371 

6613. 

14807. 

114.476 

6764. 

14807. 

W  CONV 

TOTAL 

TOTAL 

T  (Kg) 

WT  (Kg) 

COST  (M$) 

37394. 

524774. 

923. 

36419. 

270694. 

592. 

36177. 

193297. 

491. 

36100. 

155702. 

442. 

36087. 

133467. 

413. 

36108. 

118781. 

394. 

36148. 

108367. 

380. 

36202. 

100608. 

370. 

36266. 

94615. 

362. 

36338. 

89855. 

356. 

36417. 

85992. 

352. 

36502. 

82803. 

348. 

36593. 

80131. 

344. 

36689. 

77868. 

342. 

36790. 

75933. 

339. 

36897. 

74264. 

337. 

37008. 

72817. 

336. 

37124. 

71554. 

334. 

37246. 

70441. 

333. 

37372. 

69476. 

332. 

37504. 

68621. 

331. 

37642. 

67866. 

330. 

37785. 

67200. 

330. 

37933. 

66612. 

329. 

38088. 

66094. 

329. 

38249. 

65802. 

329. 

38418. 

65823. 

329. 

38249. 

65802. 

329. 

THERMODYNAMIC  PARAMETERS 
WEAPON  FLOW  RATE-Kg/s 
TURBINE  FLOW  RATE-Kg/s 
TURBINE  HYDROGEN  FLOW  RATE-kg/S 
TURBINE  OXYGEN  FLOW  RATE-kg/S 
HYDROGEN/OXYGEN  RATIO 
TURBINE  INLET  TEMP-K 
TURBINE  OUTLET  TEMP-K 
PUMP  POWER-MW 
COMBUSTOR  INLET  TEMP-K 
PRESSURE  RATIO 


10.71666813 

19.30736351 

10.62598038 

8.68138409 

1.22399616 

1150.00000000 

469.77078247 

0.88329923 

369.64312744 

114.47557831 


HYDROGEN  STORAGE  PARAMETERS 
TURBINE 
TANK  VOLUME-m3 

REFRIGERATION  SYSTEM  POWER-MW 
HYDROGEN  WEIGHT-Kg 
TANK  WEIGHT-Kg 
INSULATION  WEIGHT-kg 
REFRIG  SYST  WEIGHT-KG 
METEOROID  SHIELD  WEIGHT-kg 


107.75642395 

7.03074411E-03 

7650.70605469 

218.74552917 

351.53720093 

482.85034180 

5988.04345703 
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TABLE  5  (cont.) 


WEAPON 

TANK  VOLUME-m3 

REFRIGERATION  SYSTEM  POWER-MW 
HYDROGEN  WEIGHT-Kg 
TANK  WEIGHT-Kg 
INSULATION  WEIGHT-kg 
REFRIG  SYST  WEIGHT-KG 
METEOROID  SHIELD  WEIGHT-kg 


108.67607117 

7.07070949E-03 

7716.00097656 

220.61242676 

353.53549194 

484.93896484 

6031.96728516 


OXYGEN  PARAMETERS 
TANK  VOLUME-m3 

REFRIGERATION  SYSTEM  POWER-MW 
OXYGEN  WEIGHT-Kg 
TANK  WEIGHT-Kg 
INSULATION  WEIGHT-kg 
REFRIGERATION  SYSTEM  WEIGHT-Kg 
METEOROID  SHIELD  WEIGHT-kg 


5.47337675 

2.40834124E-04 

6250.59619141 

11.11095428 

23.92803383 

16.51741409 

461.61868286 


POWER  CONVERSION  PARAMETERS 
TURBINE  STAGES 
TURBINE  SPEED-RPM 
TURBINE  EFFICIENCY 
TURBINE  PRESSURE  RATIO 
TURBINE  WEIGHT-Kg 
GENERATOR  WEIGHT-Kg 
FLYWHEEL  WEIGHT-Kg 
POWER  CONDITIONING  WEIGHT-Kg 
NOZZLE  OUTLET  VELOCITY-m/s 


16 

10000.00000000 

0.82368428 

114.47557831 

4556.61767578 

10614.64550781 

3077.87011719 

20000.00000000 

1934.09130859 


WEIGHT  AND  COST  SUMMARY 

POWER  CONVERSION  WEIGHT-Kg 
HYDROGEN  SUBSYSTEM  WEIGHT-Kg 
OXYGEN  SUBSYSTEM  WEIGHT-Kg 
MISC.  WEIGHT-Kg 
TOTAL  WEIGHT-Kg 


38249.13281250 

14807.05468750 

6763.77148437 

5981.99609375 

65801.95312500 


POWER  CONVERSION  COST-M$ 
HYDROGEN  SUBSYSTEM  COST-M$ 
OXYGEN  SUBSYSTEM  COST-M$ 
FLYWHEEL  COST-M$ 

LAUNCH  COST-M$ 

TOTAL  COST-M$ 

Execution  terminated  :  0 


240.95716858 

0.62265307 

1.90524738E-02 

3.07787037 

84.22650146 

328.90322876 
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TABLE  6 

Combustion  Powered  Burst  System 
Weight  and  Cost  Algorithm  Summary 


Turbine  See  Appendix  C 


Generator  0.1  kg/kW  (Gerry,  1985) 


Power  Conditioning  0.2  kg/kW  (estimate) 


Flywheel  10  kg/kWh  (Bents,  1984) 


Hydrogen  Storage  Wh  “  mT 

Subsystem  V  =  W^/p 

Wt  -  2.03  V 

Wi  =  15.5  V“-667 

Wj.  -  12.6  VO  66^+109x10** 

[0. 0468 (log(2.78V**0. 667) )**2, 


$ 2000/kg 
(Gerry,  1985) 

$ 3000/kg 
(Gerry,  1985) 

$10, 000/kg 
(estimate) 

$1000/kg 

(estimate) 

hydrogen  $2 8/kg 
(Bents,  1984) 
tank  $ 20/kg 
(estimate) 
refrigeration 
9]  $640/kg  (est.) 


Ws  =  107  VO -86 

insulation 


Wh  =  hydrogen  mass,  kg 


$100/kg  (est.) 


m  =  hydrogen  flow  rate,  kg/s  shield  $10/kg 

T  =  operating  time,  s  (estimate) 

V  =  volume  of  hydrogen,  m^ 
p  -  density  of  liquid  H2,  71  kg/m^ 

Wt  =  tank  weight,  kg 
Wi  ■=  insulation  weight,  kg 
Wj.  =  refrigeration  weight,  kg 
Wg  —  meteoroid  shield  weight,  kg 
See  Appendix  D 
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TABLE  6  (cont.) 

Combustion  Powered  Burst  System 
Weight  and  Cost  Algorithm  Summary 


Component 


Weight 


Cost 


Oxygen  Storage 
Subsystem 


W^  «  mT 
Vo  -  Wo/P 
Wt  =  2.03  V„ 

Wi  -  7.7  Vo°-667 


W 


667 


.25 

■  \ 

hydrogen 
X  hydrogen  refrigerator  weight 


Same  as  for 
hydrogen  but 
oxygen  cost 
$0.2/kg 
(Gerry,  1985) 


Wg  =  107 

Wq  =  oxygen  mass,  kg 


m  oxygen  flow  rate,  kg/s 

T  =*  operating  time,  s 

Vq  =  volume  of  oxygen,  m^ 

p  =  density  of  oxygen,  1142  kg/m 

Wfc  =  tank  weight,  kg 

Wi  insulation  weight,  kg 


Oxygen  Storage 
Subsystem  cont. 


Wr  =  refrigeration  weight,  kg 
Wg  *  meteoroid  shield  weight,  kg 
See  Appendix  D 


Miscellaneous 


10%  of  subtotal 


Launch  cost 


3 


$1280/kg 
(Aviation  Week, 
1985)  and  heavy 
lift  shuttle) 


-23- 


BRAYTON  CYCLE  CONTINUOUS  POWER  SYSTEM 


The  two  previously  discussed  systems,  the  reactor  powered  and  combustion 
powered  burst  systems,  both  provided  a  high  power  level  for  a  short  time 
during  a  battle  engagement.  They  consruned  hydrogen  or  hydrogen  and 
oxygen  for  cooling  and/or  fuel.  Continuous  power  systems,  on  the  other 
hand,  must  run  for  long  periods  of  time  and  cannot  use  expendable  fuels 
and  coolants,  since  the  quantities  needed  for  long  term  operation  are 
prohibitive.  Because  of  this,  we  consider  closed  thermodynamic  cycle 
nuclear  powered,  generation  systems  to  be  the  most  likely  power  source 
for  continuous  multimegawatt  space  power.  These  systems  would  use 
thermal  radiators,  not  expendable  coolants,  to  remove  waste  heat. 

Figure  3  illustrates  a  closed  Bray ton  cycle  power  system.  A  50%  by  mass 
helium  xenon  mixture  is  used  as  the  working  fluid.  It  is  compressed, 
heated  at  constant  pressure  by  a  gas  cooled  reactor,  expanded  by  a 
turbine,  and  cooled  at  constant  pressure  by  a  radiator  to  complete  its 
cycle.  Shaft  power  from  the  turbine  drives  the  compressor  and 
generator . 

Part  of  the  flow  leaving  the  compressor  may  be  diverted  to  the  turbine 
for  blade  cooling.  Algorithms  for  calculating  the  quantity  of  blade 
coolant  required  are  given  in  Appendix  C. 

Generated  electrical  power  is  converted  by  a  power  conditioning  unit 
into  a  form  that  can  be  used  by  the  platform's  payload.  Waste  heat  from 
the  generator  and  power  conditioning  unit  is  removed  by  a  low 
temperature  radiator. 

The  efficiency  of  the  Brayton  Power  Conversion  cycle  depends  on  cycle 
temperatures  and  the  efficiencies,  of  the  turbine  and  compressor, 


_  A(Tt:i*Tt-o)  "  m(Tf’o*T(^o)  ■  (M+m)  (Tcq-Tq]^/ 
Vcyc  ~ 

M(Tti-Tco) 


(19) 


Vcyc 


1 


fl  +  s’ 

Tto  -...Tel  1 

1  M  . 

.  Tti  -  Tco  . 

(20) 


where,  M  is  the  flow  rate  through  the  turbine  excluding  blade  coolant, 
ih  is  the  blade  coolant  flow  rate, 

Tti  is  turbine  inlet  temperature, 

Tto  Is  turbine  outlet  temperature, 

Tco  Is  compressor  outlet  temperature, 

Tci  is  compressor  inlet  temperature, 
ryt  is  turbine  efficiency, 
r/c  is  compressor  efficiency,  and 
r)cyc  ^s  cycle  efficiency. 
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POWER 

CONDIT. 


FIGURE  3.  Brayton  System. 
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The  numerator  in  Equation  (19)  is  the  power  extracted  by  the  turbine 
minus  that  needed  by  the  compressor.  It  is  the  net  shaft  power 
delivered  by  the  turbine.  Notice  that  this  term  includes  the  enthalpy 
contribution  from  the  cooling  fluid  which  can  either  add  to  or  subtract 
from  turbine  power,  depending  on  the  relative  values  of  T^o  Tco- 

The  denominator  is  the  power  required  of  the  reactor.  Note  in  later 
equations  that  the  heat  generated  by  the  inefficiency  of  the  compressor 
reduces  the  energy  needed  from  the  reactor. 

The  turbine  and  compressor  outlet  temperatures  are  given  by  Equations 
(21)  and  (22). 

Tto  =  (1  -  »?t  +  »?tPr^/''''^)Tti  . 


Tco 


1 


_ 1 _ 

r?^[P^(l  +  AP)]^'^'^'^ 


(22) 


where  P^  is  the  turbine's  pressure  ratio, 

7  is  the  working  gas's  specific  heat  ratio,  and 

AP  is  the  pressure  loss  across  the  reactor  and  radiator  as  a 

fraction  of  the  turbine's  inlet  pressure. 

Because  of  the  AP  term,  the  cycle  efficiency  depends  on  pressure  losses 
in  the  system. 

In  the  cycle  analysis  T^i  is  specified,  but  Pj;  and  T^i  are  not.  They 
are  found  by  iterating  the  values  of  each  to  minimize  system  weight. 
That  is,  pressure  and  temperature  (Tti/Td)  ratios  are  optimized  to 
minimize  system  weight. 


System  efficiency  also  depends  on  generator  and  power  conditioning  unit 
efficiencies. 

(23 

^sys  ~  VcycV^V^c  > 

where  rig  is  generator  efficiency,  and 

»7pc  is  power  conditioning  efficiency. 

The  thermal  power  required  of  the  reactor  is  found  by  dividing  the 
required  electrical  power  by  system  efficiency. 

Pth  “  Pe/^sys  • 

The  power  system's  flow  rate  is  found  using  Equation  (25). 


M  —  Pth/^p  C^ti  ■  "Pco) 

M  is  the  flow  rate  and  Cp  is  the  fluid's  specific  heat. 


(25) 
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As  mentioned  above,  pressure  and  temperature  ratios  are  optimized  to  get 
minimum  system  weight.  As  the  pressure  ratio  increases,  cycle 
efficiency  increases  and  the  power  dumped  by  the  radiator  decreases; 
however,  the  radiator's  temperature  also  decreases.  These  two  things 
combine  to  give  radiator  weight  a  "U”  shaped  function  of  pressure  ratio. 
The  weights  of  other  components  enter  in  various  ways,  but  the  end 
result  is  that  there  is  a  pressure  ratio  which  minimizes  system  weight. 
Also,  as  temperature  ratio  decreases,  cycle  efficiency  decreases  and 
radiator  temperature  increases,  again  producing  a  ”U**  shaped  function  of 
temperature  ratio,  and  again  there  is  an  optimum  value  of  temperature 
ratio.  Actually,  pressure  and  temperature  ratios  have  to  be  iterated 
together  because  the  optimum  value  of  one  depends  on  the  value  of  the 
other. 

Our  Brayton  system  model  comprises  the  above  performance  algorithms,  an 
optimization  procedure,  and  algorithms  for  the  weight  of  each  component 
(discussed  later) . 

The  model  was  "computerized"  for  use  on  an  IBM-AT  personal  computer.  It 
is  written  in  FORTRAN  and  is  interactive.  Table  7  shows  a  typical  input 
sequence.  Default  values  for  all  parameters  except  power  level  and 
operation  time  are  supplied.  The  default  values  may  be  easily  changed. 
The  underlined  values  were  supplied  by  the  user.  Table  8  is  the  output. 
The  first  part  shows  the  pressure  and  temperature  ratio  optimization. 

The  second  part  shows  selected  parameter  values  for  the  optimum  system. 

Table  9  summarizes  the  weight  and  cost  algorithms  used  in  the  model. 

The  reactor  model  was  formulated  by  A1  Marshall  of  Sandia  National 
Laboratories  and  is  described  in  Marshall  (1986) .  It  uses  a  specific 
power  density  value  as  part  of  its  fuel  mass  calculation.  Power  density 
calculation  is  described  in  Appendix  A.  We  have  augmented  Marshall's 
shield  model  by  adding  a  wrap-around  shield.  This  is  described  in 
Appendix  B , 

The  turbine  model  is  described  in  Appendix  C,  and  the  compressor  weight 
algorithm  is  summarized  in  Table  9.  The  generator,  see  Table  9,  has  an 
efficiency  of  95%. 

Radiator  area  depends  on  the  quantity  of  heat  to  be  dissipated  and  on 
the  radiator's  temperature.  The  algorithms  that  calculate  radiator  area 
are  developed  in  Appendix  E,  and  radiator  weights  are  summarized  in 
Table  9. 

Power  conditioning  weight  depends  strongly  on  the  load  to  be  powered, 
but  we  have  selected  0.2  kg/kW  as  a  place  holder  until  better  estimates 
are  made.  Power  conditioning  is  assumed  to  have  a  95%  efficiency. 

The  power  conditioning  unit  and  generator  must  be  cooled  by  a  radiator. 
We  assume  that  it  is  isothermal  at  500  K  and  its  area  can  be  calculated 
like  that  of  a  Rankine  cycle  radiator. 
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TABLE  7 


THIS  PROGRAM  MODELS  A  BRAYTON  CYCLE  SPACE  POWER 
SYSTEM.  THE  SYSTEM  CONSISTS  OF  A  HELIUM  COOLED 
REACTOR,  A  TURBINE,  A  COMPRESSOR,  A  GENERATOR 
AND  A  RADIATOR.  THE  SYSTEM  OPERATING  PARAMETERS, 
TEMPERATURE  RATIO  AND  PRESSURE  RATIO,  ARE  OPTIMIZED 
TO  GET  EITHER  MINIMUM  SYSTEM  WEIGHT  OR  COST. 


ENTER  VALUES  FOR  ELECTRICAL  POWER  .IN  MW  AND 
OPERATING  TIME  IN  HOURS.  10  87600 

ENTER  "WEIGHT"  OR  "COST"  AS  PARAMETER  TO  BE 
MINIMIZED  (IN  CAPITAL  LETTERS) .  WEIGHT 

THE  FOLLOWING  ARE  DEFAULT  PARAMETER  VALUES . 
THEY  MAY  BE  CHANGED  IF  YOU  WISH. 


CYCLE  PARAMETERS : 

1.  1200.00000000 

2.  2.00000000 

3.  O.OOOOOOOOE-01 

4.  0.85000002 

5.  0.94999999 

6.  0.94999999 


TURBINE  INLET  TEMP,  K 
TURBINE  INLET  PRESSURE,  MPa 
NOT  USED 

COMPRESSOR  EFFICIENCY 

GENERATOR  EFFICIENCY 

POWER  CONDITIONING  EFFICIENCY 


REACTOR  PARAMETERS: 

7.  0.93000001 

8.  18.00000000 

9.  1.00000000 

10.  0.12200000 

11.  2.50000004E-02 

12.  2.99999993E-02 

13.  0.25000000 

14.  O.OOOOOOOOE-01 

15.  O.OOOOOOOOE-01 

16.  1.00000000 

17.  4.21999979 

18.  5560.00000000 

19.  2.70000007E-02 

20.  9.00000036E-02 

21.  4.99999992E+16 

22.  25.00000000 

23.  15.00000000 

24.  2.00000000 

25.  l.OOOOOOOOE+07 


FRACTIONAL  FUEL  ENRICHMENT 
CRITICAL  COMPACT  MASS,  Kg 
CRITICAL  MASS  CORRECTION  FACT 
FUEL  +  MODERATOR  VOL  FRACT 
FUEL  BED  LENGTH,  m 
REACTOR  PRESSUR  DROP,  FRAC  OF  PIN 
FUEL  BURNUP  FRACTION  LIMIT 
MODERATOR-TO-FUEL  RATIO 
MODERATOR  MOLECULAR  WEIGHT 
MODERATOR  DENSITY,  Kg/m3 
STRUCTURE  TO  FUEL  &  MOD  RATIO 
STRUCTURE  DENSITY,  kg/m3 
CORE  REMOVAL  X-SECTION,  cm-1 
CORE  GAMMA  ATTEN  X-SECTION,  cm-l 
ALLOWED  PAYLOAD  NEUTRON  DOSE,  nvt 
PAYLOAD  SEPARATION  DISTANCE,  10 
PROTECTION  CONE  HALF  ANGLE,  DEG 
NEUTRON  SHIELD  MATL-B4C=1,  LIH=2 
ALLOWED  PAYLOAD  GAMMA  DOSE,  R 


TURBINE  PARAMETERS: 

26.  4.00000000  NUMBER  OF  TURBINES 

27.  1.00000000  TURBINE  MATERIAL  1-SUPERALLOY, 

2 -CARBON 

28.  1200.00000000  MAXIMUM  BLADE  TEMPERATURE-K 
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TABLE  7  (cont.) 


29.  10000.00000000  TURBINE  SPEED-RPM 

30.  2.00000000  WORK  COEFFICIENT 


RADIATOR  PARAMETERS: 

31.  0.88000000  RADIATOR  EMITTANCE 


ENTER  THE  NUMBER  OF  PARAMETERS  YOU  WISH  TO  CHANGE. 
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TABLE  8 


PRESSURE 

TEMP 

CYCLE 

RCT+SHLD 

POW  CONV 

RADIATOR 

TOTAL 

TOTAL 

RATIO 

RATIO 

EFFIC 

WT  (Kg) 

WT  (Kg) 

WT  (Kg) 

WT  (Kg) 

COST  (M$) 

1.800 

2.527 

0.137 

18872. 

17810. 

108129. 

159291. 

357. 

2.600 

2.924 

0.222 

12786. 

14654. 

101952. 

142332. 

301. 

3.400 

3.253 

0.280 

11451. 

13982. 

115871. 

155434. 

313. 

2.200 

2.736 

0.184 

14849. 

15907. 

100438. 

144313. 

315. 

3.000 

3.095 

0.254 

11744. 

13737. 

107592. 

146380. 

301. 

2.400 

2.833 

0.204 

13677. 

15226. 

100508. 

142352. 

306. 

2.800 

3.012 

0.239 

12083. 

14164. 

104386. 

143696. 

299. 

2.500 

2.879 

0.213 

13203. 

14929. 

101087. 

142141. 

303. 

2.700 

2.968 

0.231 

12415. 

14400. 

103062. 

142865. 

300. 

2.500 

1.727 

-0.246 

2.500 

2.015 

0.078 

47493. 

106478. 

160475. 

345891. 

910. 

2.500 

2.303 

0.157 

19759. 

28959. 

89774. 

152341. 

371. 

2.500 

2.591 

0.193 

14297. 

18700. 

89085. 

134290. 

303. 

2.500 

2.879 

0.213 

13203. 

14929. 

101087. 

142141. 

303. 

2.500 

2.447 

0.178 

16215. 

22318. 

87074. 

138167. 

323. 

2.500 

2.735 

0.204 

13660. 

16451. 

93976. 

136495. 

299. 

2.500 

2.519 

0.186 

15012. 

20272. 

87648. 

135225. 

310. 

2.500 

2.663 

0.199 

13949. 

17457. 

91227. 

134896. 

300. 

2.500 

2.555 

0.190 

14501. 

19438. 

88270. 

134429. 

305. 

2.500 

2.627 

0.196 

14115. 

18044. 

90075. 

134457. 

301. 

1.800 

2.274 

0.121 

20833. 

23136. 

101157. 

159639. 

375. 

2.600 

2.632 

0.201 

13822. 

18294. 

89256. 

133509. 

299. 

3.400 

2.928 

0.257 

11705. 

16656. 

96772. 

137646. 

293. 

2.200 

2.463 

0.165 

16175. 

20168. 

90468. 

139493. 

321. 

3.000 

2.786 

0.231 

12380. 

16952. 

91924. 

133382. 

290. 

2.800 

2.711 

0.217 

13025. 

17574. 

90250. 

132934. 

294. 

3.200 

2.858 

0.245 

11854. 

16795. 

94130. 

135057. 

290. 

2.700 

2.672 

0.209 

13401. 

17920. 

89655. 

133074. 

296. 

2.900 

2.749 

0.224 

12686. 

17252. 

91013. 

133047. 

292. 

2.800 

1.807 

-0.269 

2.800 

2.108 

0.094 

41434. 

91382. 

147747. 

308620. 

803. 

2.800 

2.409 

0.179 

17971. 

26717. 

88608. 

146626. 

351. 

2.800 

2.711 

0.217 

13025. 

17^74. 

90250. 

132934. 

294. 

2.800 

3.012 

0.239 

12083. 

14164. 

104386. 

143696. 

299. 

2.800 

2.560 

0.201 

14828. 

20819. 

87209. 

135142. 

310. 

2.800 

2.861 

0.229 

12477. 

15544. 

96150. 

136588. 

293. 

2.800 

2.635 

0.210 

13754. 

18988. 

88311. 

133157. 

300. 

2.800 

2.786 

0.224 

12726. 

16453. 

92889. 

134276. 

292. 

2.800 

2.673 

0.214 

13296. 

18238. 

89186. 

132792. 

296. 

2.800 

2.748 

0.220 

12869. 

16983. 

91488. 

133473. 

293. 

2.800 

2.673 

0.214 

13296. 

18238. 

89186. 

132792. 

296. 

CYCLE  PARAMETERS 

CYCLE  EFFICIENCY 
THERMAL  POWER-MW 
MASS  FLOW  RATE-Kg/S 
TURBINE  INLET  TEMP-K 
TURBINE  OUTLET  TEMP-K 
COMPRESSOR  INLET  TEMP-K 
COMPRESSOR  OUTLET  TEMP-K 
PRESSURE  RATIO 


0.21353653 

51.88963699 

41.14289474 

1200.00000000 

820.02722168 

448.94909668 

728.16864014 

2.79999995 
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TABLE  8  (cont.) 


REACTOR  PARAMETERS 
BURNUP  MASS-Kg 
INITIAL  CRITICAL  MASS-Kg 
END  CRITICAL  MASS-Kg 
TOTAL  BRNUP+CRIT  MASS-Kg 
SPECIFIC  POWER-W/Kg 
MASS  FOR  SPECIFIC  POWER  LIM-Kg 
MASS  FOR  ALLOWED  BURNUP-Kg 
FUEL  MASS-Kg 
MODERATOR  MASS -Kg 
STRUCTURE  MASS-Kg 
REFLECTOR  MASS-Kg 
PRESSURE  VESSEL  MASS-Kg 
MISCELLANEOUS  MASS-Kg 
TOTAL  REACTOR  MASS-Kg 
NEUTRON  SHIELD  THICKNESS-Itl 
NEUTRON  SHIELD  MASS-Kg 
GAMMA  SHIELD  THICKNESS-in 
GAMMA  SHIELD  MASS-Kg 
TOTAL  SHIELD  MASS-Kg 


RADIATOR  PARAMETERS 
INLET  TEMPERATURE-K 
OUTLET  TEMPERATURE-K 
HIGH  TEMPERATURE  AREA-ltl2 
MEDIUM  TEMPERATURE  AREA-in2 
LOW  TEMPERATURE  AREA-m2 
TOTAL  AREA-m2 
TOTAL  WEIGHT-Kg 


POWER  CONVERSION  PARAMETERS 
TURBINE  STAGES 
TURBINE  SPEED-RPM 
TURBINE  EFFICIENCY 
TURBINE  PRESSURE  RATIO 
TURBINE  COOLANT  FLOW  RATE-kg/s 
TURBINE  WEIGHT-Kg 
COMPRESSOR  WEIGHT-Kg 
GENERATOR  WEIGHT-Kg 
POWER  CONDITIONING  WEIGHT-Kg 
GEN  &  PC  RADIATOR  WEIGHT-Kg 
GEN  &  PC  RADIATOR  AREA-in2 
TOTAL  WEIGHT-Kg 


=  212.02217102 
=  697.88366699 
=  712.42236328 
=  924.44451904 
=  65079.46484375 
=  1036.52551270 
=  1025.33923340 
=  1036.52551270 
=  O.OOOOOPOOE-01 
=  2268.67578125 
=  1144.35083008 
=  339.70037842 
=  1036.52551270 
=  5825.77783203 
=  0.22490446 
=  877.93627930 
=  7.39577487E-02 
=  6592.65576172 
=  7470.59179687 


770.83776855 

426.50921631 

O.OOOOOOOOE-01 

1210.29101562 

7576.68359375 

8786.97460937 

89185.77343750 


=  14 

=  8266.25000000 

=  0.93672222 

=  2.79999995 

=  O.OOOOOOOOE-01 

=  3221.94262695 

=  9665.82812500 

=  1052.63159180 

=  2000.00000000 

=  2297.35815430 

=  367.57730103 

=  18237.75976562 
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TABLE  8  (cont.) 


WEIGHT  AND  COST  SUMMARY 
REACTOR  WEIGHT-Kg 
SHIELD  WEIGHT-Kg 
POWER  CONVERSION  WEIGHT-Kg 
RADIATOR  WEIGHT-Kg 
MISCELLANEOUS  WEIGHT-Kg 
TOTAL  WEIGHT-Kg 


=  5825.77783203 
=  7470.59179687 
=  18237.75976562 
=  89185.77343750 
=  12071.99023437 
=  132791.89062500 


REACTOR+SHIELD  COST-M$ 
POWER  CONVERSION  COST-M$ 
RADIATOR  COST-M$ 

LAUNCH  COST-M$ 

TOTAL  COST-M$ 

Execution  terminated  :  0 


58.50402451 

49.39291000 

17.83715439 

.169.97361755 

295.70770264 
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TABLE  9 

Bray ton  System 

Weight  and  Cost  Algorithm  Summary 


Component _ _ Weight _ _ Cost 


Reactor  and  Shield 

Marshall's  algorithm 
(Marshall,  1986) 

$4400/kg 

(estimate) 

Turbine 

See  Appendix  C 

$ 2000/kg 
(Gerry,  1985) 

Generator 

0.1  kg/kW  (Gerry,  1985) 

$3000/kg 
(Gerry,  1985) 

Power  Conditioning 

0.2  kg/kW  (estimate) 

$10,000Ag 

(estimate) 

Radiator 

12  kg/m2  for  temp  >  1000  K 

8  kg/m^  for  650  <  temp  <  1000  K 

5  kg/m^  for  temp  <  650  K 
(NASA  estimates) 

$ 200/kg 
(estimate) 

multiply  by  1.25  for  meteoroid  loss 
(estimate) 

multiply  by  1.5  for  heat  exchanger 
(estimate) 

See  Appendix  F  for  area  calculation 

Compressor 

3  times  turbine  (estimate) 

$2000/kg 
(Gerry,  1985) 

PC  6c  Generator 

Radiator 

5  kg/m2 

(NASA  estimate) 

$200/kg 

(estimate) 

add  25%  for  meteoroid  loss 
(estimate) 

Miscellaneous 

10%  of  subtotal 

Launch  cost 

$1280/kg 
(Aviation  Week, 
1985,  and  heavy 
lift  shuttle) 
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RANKTHK  CYCLE  CONTINUOUS  POWER  SYSTEM 


Like  the  nuclear  reactor  powered  Brayton  cycle  system,  the  nuclear 
reactor  powered  liquid  metal  Rankine  cycle  system  can  potentially 
provide  continuous  power  for  several  years .  Because  of  its  long 
operation  time,  expendables  cannot  be  practically  used.  Cycle  waste 
heat  is  rejected  by  a  space  radiator,  and  this  is  where  the  Rankine 
system  has  a  significant  advantage  over  the  Brayton  system.  The  Rankine 
system  rejects  heat  from  a  condensing  working  fluid;  thus,  its  radiator 
is  nearly  isothermal  and  radiates  much  more  heat  per  unit  area  than  a 
comparable  Brayton  radiator.  Since  the  Brayton  cycle  rejects  sensible 
heat  from  its  working  fluid,  its  radiator  experiences  a  substantial 
temperature  drop  from  the  inlet  end  to  the  outlet  end.  Thus,  Rankine 
radiators  use  much  less  area  than  Brayton  radiators  to  reject  the  same 
quantity  of  heat.  A  major  disadvantage  of  Rankine  systems  is  that  the 
technology  associated  with  using  two  phase  alkaline  metals  as  working 
fluids  in  space  is  not  well  developed. 

We  have  modeled  two  types  of  Rankine  systems,  direct  and  indirect,  shown 
in  Figures  4  and  5  respectively.  The  direct  system  boils  potassium  in  a 
liquid  metal  cooled  reactor  and  sends  the  saturated  vapor  to  a  turbine 
for  power  generation.  Since  the  fluid  leaving  the  "hot"  end  of  the 
reactor  is  unlikely  to  be  pure  vapor,  a  separator  is  used  to  separate 
the  saturated  vapor  from  its  accompanying  liquid.  The  liquid  is 
recirculated  to  the  "cold"  end  of  the  reactor. 

The  indirect  Rankine  system  uses  liquid  lithium  as  a  reactor  coolant. 

It  transfers  heat,  via  heat  exchangers,  to  the  potassium  working  fluid. 
Our  model  of  this  system  allows  for  nearly  any  level  of  superheat  as 
indicated  by  the  temperature  entropy  diagrams  in  Figures  6a,  6b,  and  6c. 

Our  model  for  the  Rankine  cycle  is  an  approximate  one,  but  it  gives  very 
accurate  results.  We  divide  the  ideal  cycle  into  three  zones  -- 
preheat,  boil,  and  superheat  --as  shown  in  Figure  7.  Then  we  calculate 
the  Carnot  efficiency  for  each  zone  based  on  the  average  temperature  in 
that  zone.  Finally,  we  weight  the  three  efficiencies  by  the 
corresponding  enthalpy  added  to  the  working  fluid  in  that  zone  to  get 
the  ideal  cycle  efficiency,  and  we  multiply  the  ideal  cycle  efficiency 
by  turbine  efficiency  to  get  the  cycle  efficiency.  We  have  compared 
this  method  to  one  using  real  enthalpy  values  for  several  cases  and  find 
that  our  approximation  is  usually  within  2  to  3  percent  of  the 
rigorously  derived  values.  By  2  to  3  percent,  we  mean  that  if  the 
rigorous  value  is  20%,  the  approximate  value  will  be  19.4%,  not  17%. 

The  following  steps  will  show  this  method  in  more  detail. 
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PUMP 


FIGURE  4.  Direct  Rankine  System 


POWER 

CONDITION 

ING 
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PUMP 


FIGURE  5 .  Indirect  Rankine  System 


POWER 

CONDIT. 
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(a) 


FIGURE  6a.  Direct  Rankine  Cycle  and 
Indirect  Cycle  Without  Superheat 


FIGURE  6b.  Indirect  Rankine  Cycle 
With  Some  Superheat 


FIGURE  6c.  Indirect  Rankine  Cycle 
With  Superheat 


ENTROPY 
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FIGURE  7.  Ideal  Rankine  Cycle 
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1.  Preheat  zone  --  The  equations  for  this  zone  are  as  follows: 

.5(Tb  +  TJ  -  Te 

r?p  -  -  ,  (26) 

.5(Tb  +  TJ 

where  t/p  is  the  efficiency  in  the  preheat  zone, 

Tb  is  the  boiling  temperature,  and 
Tjj  Is  the  condensing  temperature. 

The  enthalpy  added  in  this  zone  is  given  in  equation  27 : 

hp  =  Cpi(Tb  -  T,)  ,  (27) 

where  Cp^^  is  the  potassium  liquid's  specific  heat  (830  J/kg)  . 

2.  Boiling  zone  --  The  equations  for  this  zone  are  as  follows: 

Tb  -  Te 

- ,  (28) 

Tb 

where  is  the  efficiency  in  the  preheat  zone. 

The  boiling  enthalpy  depends  on  boiling  temperature.  Fitting  potassium 
data  from  Reynolds  (1979)  gives  the  following  result: 

hfg  -  2.57  X  106  .  640Tb  »  (29) 

where  hfg  is  the  enthalpy  of  vaporization  for  potassium. 


3.  Superheat  zone  --  The  equations  for  this  zone  are  as  follows: 

.5(T3  +  Tb)  -  Te 

Vs  =  -  ,  (30) 

.5(T3  +  Tb) 

where  is  the  efficiency  in  the  superheat  zone,  and 

Tg  is  the  peak,  or  turbine  inlet  temperature. 


The  superheat  enthalpy  is  given  by: 

K  -  Cpg(T3  -  Tb)  ,  (31) 

where  hg  is  the  superheat  enthalpy,  and 
Cpg  is  the  gas  specific  heat. 
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Cpg  was  also  found  by  fitting  data  in  Reynolds  (1979). 

Cpg  =  1.7Tx  -  .05Ts  -  .0005TxTs  .  (^2) 

Tx  is  the  saturation  temperature  of  the  gas .  That  is ,  it  is  the 
temperature  where  the  gas  will  become  a  saturated  vapor  if  it  is 
expanded  isentropically .  It  was  also  found  by  fitting  data  from 
Reynolds  (1979). 

Tx  ==  295  e2*6/S  ^  (33) 


S  =  Sq  +  ln(Tg/T|3) /1 . 06 


So  =  2.6/[ln(Tb)  -  5.687]  .  (35) 

In  these  equations,  S  and  S©  are  the  entropy  at  Tg  and  the  entropy  at  Tb 
respectively. 


Equations  (26)  through  (35)  and  turbine  efficiency  (see  Appendix  C)  are 
all  we  need  to  estimate  cycle  efficiency: 


where 


Vcyc  ■“ 


hp»?p  +  hb»/b  + 
nt - - 

hp  +  hb  +  hg 


ricyc  is  cycle  efficiency,  and 
r?t  is  turbine  efficiency. 


(36) 


For  an  example,'  let  Tg  =  1200  K,  Tb  =  1109  K,  Tc  =  900  K,  and 
We  get  the  following: 


Tx  =  1057, 

hg  =  100.4  kJ/kg, 

*?s  =  -220, 

hb  -  I860  kJ/kg, 

rih  =  -1885, 

hp  -  173.5  kJ/kg, 

»;p  =  .  104 ,  and 
’Jcyc  “  16.48%. 

Using  enthalpy  and  entropy  values  from  Reynolds  (1979)  and  Figure  7: 


hi  =  79.3  kJ/kg, 

h3  =  2225.1  kJ/kg,  and 

hs  -  1087.8  kJ/kg. 

S3  =  S4  -  20418  kJ/kg  K, 
Si  =  .0933  kJ/kg  K,  and 
Ss  =  2.3250  kJ/kg  K. 
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X  =  quality  at  point  4  -  -  .8731,  and 

i>5  -  5>i 

h4  =  hi  +  x(h5-hi)  ”  1832.9  kJ/kg. 

We  also  need  pump  enthalpy,  hp,  which  is  equal  to  pressure  change  times 
specific  volume  divided  by  efficiency  (.7). 

Pi  -  .025  X  106  Pa, 

P2  =  . 2  X  10®  Pa ,  and 
Vi  =  .001437  m3/kg. 

(P2  -  Pi)vi 

hpump - 359  J/kg 

.7 


•9(h3  -  h^)  -  hpu^ 

^  eye  “16.44% 

hs  *  hi  -  hpump 

Thus,  our  approximate  model  calculates  a  cycle  efficiency  of  16.48%,  and 
a  rigorous  analysis  calculates  an  efficiency  of  16.44%  for  our  example. 

In  our  model,  we  could  specify  Tg,  which  is  both  the  superheat 
temperature  and  the  turbine  inlet  temperature,  and  the  boiler 
temperature,  Tj,,  but  we  do  not.  Instead,  we  specify  Tg  and  the  ratio  of 
turbine  inlet  pressure  to  saturation  pressure  at  the  turbine  inlet 
temperature,  and  we  use  this  ratio  to  calculate  the  boiler  temperature. 
This  is  done  using  a  derivative  of  the  Claus ius-Clapeyron  equation  which 
gives  the  relation  between  boiling  temperature  and  boiling  pressure. 

For  potassium,  this  equation  is  as  follows: 


-9760/Tboii 

Pboii  =  1320  e 


(37) 


Thus: 

Pti  1320  e 
1320 

and 

T, _ 9760  T„ 

9760  -  Tgln(Pti/Pg) 


(38) 


(39) 
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where 


Pti  is  turbine  inlet  pressure,  also  the  boiler  pressure  (Pa) 
and 

Pg  is  the  saturated  vapor  or  boiling  pressure  that  corresponds 
to  a  temperature  of  Tg. 

If  p^j^  —  Pg,  then  the  cycle  does  not  use  any  superheat.  As  Pti/l*s  gi^ows 
less  than  1,  the  quantity  of  superheat  grows. 

Turbine  inlet  and  condenser  pressures  can  be  calculated  using 
equation  37 : 


Pti 


-  1320  e 


-9760/Tb 


(40) 


Pc 


-  1320  e 


-9760/Tc 


(41) 


The  heat  which  must  be  supplied  by  the  system* s  reactor,  Qrctr>  given 
by  equation  42: 


Qrctr  Pe/  C9cyc9g9pc) 


(42) 


where  Pe  is  the  electrical  power  required  (W)  , 
fjg  is  generator  efficiency,  and 
»;pe  is  power  conditioning  efficiency. 

The  heat  rejected  by  the  radiator,  Qrad.  is  given  by  equation  43: 
Qrad  “  Qrctr  “  Pe/(9g9po) 


Turbine  power,  Pturb»  is  given  by  equation  44: 
Pturb  ~  Pe/ (9g9pc) 


(44) 


Working  fluid  flow  rate,  m,  is  given  by  equation  45: 


m  ”  Qrctr/ (^)  +  hb  +  hg) 


(45) 


This  is  also  the  reactor  flow  rate  for  the  direct  cycle,  but  hg  is  zero. 
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To  find  the  reactor  flow  rate,  for  an  indirect  cycle,  we  have  to 

make  some  assumptions  about  heat  exchanger  effectiveness.  Figure  8 
shows  a  temperature -enthalpy  plot  for  the  indirect  cycle.  Working  fluid 
enters  the  preheater  and  progresses  to  the  superheater.  Reactor  fluid 
flows  in  the  opposite  direction,  and  it  must  satisfy  "pinch  point"  • 
requirements.  That  is,  Tpp  must  be  greater  than  Tb.  Reactor  inlet  and 
outlet  temperatures  and  flow  rate  can  be  found  as  follows : 

Tb  ■  Tc 

Tpp  -  Tc  +  - ^ ,  (46) 

'Jpx 


Ts  -  Tb 

Tro  Tb  +  , 

^sx 

where  r^p^  is  preheater  effectiveness  (.9),  and 
t7sx  is  superheater  effectiveness  (.8). 

But  we  also  require  that  Tro  -  Tpp  +  1. 

.  _  ih(hb  +  hs) 

CpLidro  -  Tpp)  ’ 

where  CpLi  is  the  lithium  reactor  coolant's  specific  heat. 


(47) 


(48) 


We  also  require  that  must  not  exceed  5  kg/s  per  MW^h  generated  in  the 
reactor.  This  5  kg/s  per  MW^h  number  represents  one  of  the  highest  flow 
rates  we  have  seen  for  a  liquid  metal  cooled  reactor.  Higher  flow  rates 
may  be  practical,  but  we  have  not  yet  studied  their  implications.  If  we 
use  the  flow  rate  limit,  then  we  must  adjust  Tj-q- 


If  ihj-  >  5  kg/MW-s,  then  ihj-  =  5  kg/MW-s,  and 

Tro  ^  *^PP  ^  ^(^b  ^s)/^r^PLi  • 

The  reactor's  inlet  temperature  is  given  by  Equation  (50): 
Tri  Tpp  -  riil^p/^rCpLi  i 


(49) 

(50) 


and  reactor  pressure  is  assumed  to  be  the  saturated  liquid  pressure  for 
lithium  at  Tro' 


^react  ^  7200 


^.18100/Tro 


(51) 
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This  is  the  Claus ius-Clapeyron  derivative  equation  for  lithium. 

Condenser  temperature  is  a  variable  which  can  be  optimized.  As  it 
decreases,  the  cycle  becomes  more  efficient,  reactor  mass  decreases,  and 
the  waste  heat  which  must  be  removed  by  the  radiator  decreases.  On  the 
other  hand,  as  condenser  temperature  decreases,  radiator  temperature 
decreases  and  the  radiator  can  remove  less  heat  per  unit  area.  Thus, 
there  is  an  optimum  condenser  temperature  which  minimizes  system  mass. 

Our  Rankine  system  model  was  "computerized"  for  use  on  the  IBM-AT 
personal  computer.  It  is  written  in  FORTRAN/77  and  is  interactive.  The 
program  uses  the  above  algorithms  to  determine  system  performance  and 
parameter  values.  Then,  it  estimates  mass  and  cost  values  for  each 
component  and  for  the  system  (discussed  later) .  Table  10  shows  a 
typical  input  sequence  for  running  the  computer  program  which  we  call 
"RNKCYC."  First,  the  user  is  asked  to  specify  a  power  level  in  MWe  and 
an  operation  time  in  hours  (user  input  is  underlined  in  Table  10). 

Next,  default  parameter  values  are  listed  and  changes  are  solicited. 

When  all  of  the  changes  have  been  made,  the  program  runs  and  prints  the 
information  shown  in  Table  11. 

The  first  part  of  Table  11  shows  the  optimization  of  condenser 
temperature.  In  this  example,  the  minimum  mass  system  is  achieved  for  a 
675  K  condenser  temperature.  Notice  that  the  system  mass  values 
decrease  until  condenser  temperature  is  825  K.  Then,  mass  starts 
increasing  until  it  takes  a  sharp  drop  at  675  K  and  starts  increasing 
again.  The  drop  was  caused  by  switching  radiator  materials  from 
titanium  to  altiminiim.  The  switch  was  made  because  aluminum  is  used  when 
radiator  temperature  is  below  650  K.  (There  is  a  30  K  difference 
between  condenser  and  radiator  temperature.) 

The  second  part  of  Table  11  lists  parameter  values  for  the  optimum 
system. 

Table  12  summarizes  component  mass  algorithms  which  are  also  described 
below. 

The  model  assumes  the  use  of  a  lithium  (indirect  cycle)  or  potassium 
(direct  cycle)  cooled  reactor.  The  reactor  model  was  formulated  by  Al 
Marshall  of  Sandia  National  Laboratories  and  is  described  in  Marshall 
(1986).  As  with  other  reactors  already  discussed,  his  model  calculates 
three  fuel  mass  values  --  end-of-life  criticality,  burnup  limit,  and  the 
specific  power  fuel  requirement  --  and  uses  the  largest  of  the  three. 

To  the  fuel  mass  is  added  moderator  mass  (if  any),  structure  mass, 
reflector  mass,  pressure  vessel  mass,  a  miscellaneous  mass,  and  gamma 
and  neutron  shield  mass. 

Radiator  mass  is  calculated  using  the  mass  per  unit  area  values  in 
Table  12,  and  the  area  algorithms  in  Appendix  F.  To  the  radiator  is 
added  a  heat  exchanger  (condenser)  mass  calculated  using  the  algorithms 
in  Table  12. 
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Turbine  mass  algorithms  are  given  in  Appendix  C.  They  were  derived  from 
algorithms  developed  by  Steve  Hudson  of  Sandia  National  Laboratories . 

The  direct  cycle  uses  a  vapor  separator  and  its  mass  is  calculated  using 
the  same  algorithms  as  for  a  condenser.  The  indirect  cycle  does  not  use 
a  separator  but  it  does  use  a  heat  exchanger.  Heat  exchanger  mass 
algorithms  are  given  in  Table  12.  Alternator  mass  is  assumed  to  be 
0.1  kgAWe  and  is  based  on  a  near  term,  iron  core  type  of  alternator. 

It  has  an  efficiency  of  95%. 

Power  conditioning  mass  is  assumed  to  be  0.2  kg/kW.  This  is  a  place 
holder  and  the  real  value  will  depend  on  load  characteristics.  We 
assume  power  conditioning  is  95%  efficient.  The  alternator  and  power 
conditioning  must  be  cooled  and  we  assxame  the  use  of  a  radiator  for 
this.  It  rejects  heat  at  500  K  and  is  assumed  to  be  isothermal.  We 
multiply  its  area  (calculated  using  the  Rankine  radiator  algorithm  in 
Appendix  F)  by  1.25  to  account  for  meteoroid  losses. 

A  miscellaneous  mass  of  10%  is  added  to  account  for  structure,  piping, 
etc. 

Costs  are  estimated,  but  they  are  at  present  very  crude  and  little 
confidence  should  be  placed  in  them. 
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TABLE  10 


THIS  PROGRAM  MODELS  A  RANKINE  CYCLE  SPACE  POWER 
SYSTEM,  THE  SYSTEM  CONSISTS  OF  A  LIQUID  METAL 
COOLED  REACTOR,  A  TURBINE,  A  GENERATOR 
A  RADIATOR,  AND  POWER  CONDITIONING. 

THREE  TYPES  OF  CYCLES  ARE  COVERED: 

1.  A  DIRECT  CYCLE  WITH  NO  SUPERHEAT,  A  VAPOR 
SEPARATOR,  AND  LITHIUM  AS  A  WORKING  FLUID. 

2.  AND  AN  INDIRECT  CYCLE  WITH  SUPERHEAT,  A  HEAT- 
EXCHANGER,  POTASSIUM  AS  A  WORKING  FLUID  AND 
LITHIUM  AS  A  REACTOR  COOLANT. 


CONDENSER  TEMPERATURE  IS  OPTIMIZED 

TO  GET  EITHER  MINIMUM  SYSTEM  WEIGHT  OR  COST. 


ENTER  VALUES  FOR  ELECTRICAL  POWER  IN  MW  AND 
OPERATING  TIME  IN  HOURS.  10  87600 

ENTER  "WEIGHT”  OR  "COST”  AS  PARAMETER  TO  BE 
MINIMIZED  (IN  CAPITAL  LETTERS) .  WEIGHT 

THE  FOLLOWING  ARE  DEFAULT  PARAMETER  VALUES. 
THEY  MAY  BE  CHANGED  IF  YOU  WISH. 


CYCLE  PARAMETERS: 

1.  1200.00000000 
2.  1.00000000 

3.  1.00000000 

4.  0.94999999 

5.  0.94999999 


TURBINE  INLET  TEMP,  K 

CYCLE  TYPE  1-DIRECT,  2 -INDIRECT 

TURBINE  INLET/SATURATION  PRESSURE 

GENERATOR  EFFICIENCY 

POWER  CONDITIONING  EFFICIENCY 


REACTOR 

6. 

7. 

8. 
9. 

10. 

11. 

12. 

13. 

14. 

15. 

16. 

17. 

18. 

19. 

20. 
21. 
22. 

23. 

24. 


PARAMETERS: 

0.93000001 

18.00000000 

12170-00000000 

1.00000000 

0.28000000 

150000.00000000 

7.00000003E-02 

O.OOOOOOOOE-01 

O.OOOOOOOOE-01 

1.00000000 

1.42999995 

8300.00000000 

0.83999997 

0.23300000 

4.99999992E+16 

25.00000000 

15.00000000 

2.00000000 

l.OOOOOOOOE+07 


FRACTIONAL  FUEL  ENRICHMENT 
CRITICAL  COMPACT  MASS,  Kg 
FUEL  DENSITY,  kg/ia3 
CRITICAL  MASS  CORRECTION  FACT 
FUEL  +  MODERATOR  VOL  FRACT 
FUEL  POWER  DENSITY,  W/kg 
FUEL  BURNUP  FRACTION  LIMIT 
MODERATOR-TO-FUEL  RATIO 
MODERATOR  MOLECULAR  WEIGHT 
MODERATOR  DENSITY,  Kg/in3 
STRUCTURE  TO  FUEL  &  MOD  RATIO 
STRUCTURE  DENSITY,  kg/in3 
CORE  REMOVAL  X-SECTION,  cm-l 
CORE  GAMMA  ATTEN  X-SECT,  cm-l 
ALLOWED  PAYLOAD  NEUTRON  DOSE,  nvt 
PAYLOAD  SEPARATION  DISTANCE,  m 
PROTECTION  CONE  HALF  ANGLE,  DEG 
NEUTRON  SHIELD  MATL-B4C=1,  LIH=2 
ALLOWED  PAYLOAD  GAMMA  DOSE,  R 
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TABLE  10  (cont.) 


TURBINE  PARAMETERS 

25.  4.00000000 

26.  1.00000000 

27. '  10000.00000000 

28.  1.50000000 

29.  O.OOOOOOOOE-01 


NUMBER  OF  TURBINES 
TURBINE  MATERIAL  1-Ni,  2-Ta 
TURBINE  SPEED,  RPM 
TURBINE  WORK  COEFFICIENT 
DISK  COOL  PARAM,  O-NO,  l-YES 


RADIATOR 

30. 


PARAMETERS : 

0.88000000 


RADIATOR  EMITTANCE 


ENTER  THE  NUMBER  OF  PARAMETERS  YOU  WISH  TO  CHANGE. 


0 


TABLE  11 


TURBINE 

CONDENSER 

CYCLE 

RCT+SHLD 

TEMP  (K) 

TEMP  (K) 

EFFIC 

WT  (Kg) 

1200.000 

1125.000 

0.055 

66009. 

1200.000 

1100.000 

0.072 

51471. 

1200.000 

1075.000 

0.090 

42548. 

1200.000 

1050.000 

0.106 

36501. 

1200.000 

1025.000 

0.123 

32121. 

1200.000 

1000.000 

0.139 

28796. 

1200.000 

975.000 

0.154 

26181. 

1200.000 

950.000 

0.170 

24068. 

1200.000 

925.000 

0.185 

22323. 

1200.000 

900.000 

0.200 

20857. 

1200.000 

875.000 

0.214 

19606. 

1200.000 

850.000 

0.228 

18525. 

1200.000 

825.000 

0.242 

17582. 

1200.000 

800.000 

0.256 

16750. 

1200.000 

775.000 

0.269 

16011. 

1200.000 

750.000 

0.283 

15350. 

1200.000 

725.000 

0.296 

14754. 

1200.000 

700.000 

0.309 

14214. 

1200.000 

675.000 

0.321 

13722. 

1200.000 

650.000 

0.334 

13272. 

1200.000 

625.000 

0.346 

12858. 

1200.000 

600.000 

0.358 

12476. 

1200.000 

575.000 

0.370 

12121. 

1200.000 

550.000 

0.382 

11792. 

1200.000 

525.000 

0.394 

11484. 

1200.000 

500.000 

0.405 

11196. 

1200.000 

475.000 

0.417 

10926. 

1200.000 

450.000 

0.428 

10672. 

1200.000 

425.000 

0.440 

10431. 

1200.000 

400.000 

0.451 

10204. 

1200.000 

375.000 

0.462 

9988. 

1200.000 

350.000 

0.473 

9783. 

1200.000 

325.000 

0.485 

9587. 

1200.000 

300.000 

0.496 

9400. 

1200.000 

675.000 

0.321 

13722. 

>OW  CONV 

RADIATOR 

TOTAL 

TOTAL 

WT  (Kg) 

WT  (Kg) 

WT  (Kg) 

COST  (M$) 

25568. 

39706. 

144411. 

547. 

20769. 

32438. 

115145. 

435. 

17918. 

28326. 

97672. 

367. 

16045. 

25824. 

86208. 

321. 

14733. 

16182. 

69340. 

276. 

13774. 

15564. 

63948. 

252. 

13054. 

15248. 

59931. 

234. 

12504. 

15166. 

56912. 

220. 

12080. 

15281. 

54653. 

208. 

11755. 

15573. 

53003. 

199. 

11599. 

16034. 

51962. 

192. 

11528. 

16665. 

51390. 

187. 

11538. 

17474. 

51253. 

182. 

11628. 

18477. 

51540. 

180. 

11799. 

19698. 

52260. 

178. 

12116. 

21170. 

53499. 

177. 

12524. 

22936. 

55235. 

178. 

13036. 

25053. 

57534. 

180. 

13674. 

17248. 

49109. 

167. 

14463. 

19165. 

51591. 

170. 

15438. 

21493. 

54767. 

175. 

16643. 

24341. 

58805. 

181. 

18136. 

27862. 

63932. 

190. 

19997. 

32269. 

70464. 

201. 

22332. 

37867. 

78852. 

217. 

25285. 

45108. 

89748. 

237. 

29056. 

54683. 

104132. 

263. 

33926. 

67707. 

123535. 

299. 

40298. 

86073. 

150483. 

.  349. 

48760. 

113283. 

189472. 

421. 

60191. 

156547. 

249398. 

528. 

75935. 

233420. 

351051. 

704. 

98114. 

400527. 

559051j^ 

1047, 

130189. 

997690. 

1251007. 

2115, 

13674. 

17248. 

49109. 

167. 

CYCLE  PARAMETERS 
CYCLE  EFFICIENCY 
THERMAL  POWER-MW 
MASS  FLOW  RATE-Kg/s 
TURBINE  INLET  TEMP-K 
CONDENSER  TEMP-K 
BOILER  TEMP-K 
BOILER  PRESSURE-MPa 
CONDENSER  PRESSURE-MPa 
REACTOR  OUTLET  TEMP-K 
REACTOR  INLET  TEMP-K 
REACTOR  PRESSURE-MPa 
REACTOR  FLOW  RATE-kg/s 


0.32119268 

34.49746323 

15.41613865 

1200.00000000 

675.00000000 

1200.00000000 

0.38753614 

6.93421869E-04 

1200.00000000 

675.00000000 

0.38753614 

15.41613865 
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TABLE  11  (cont.) 


REACTOR  PARAMETERS 
BURNUP  MASS«Kg 
INITIAL  CRITICAL  MASS-Kg 
END  CRITICAL  MASS-Kg 
TOTAL  BRNUP+CRIT  MASS-Kg 
SPECIFIC  POWER-W/Kg 
MASS  FOR  SPECIFIC  POWER  LIM-Kg 
MASS  FOR  ALLOWED  BURNUP-Kg 
FUEL  MASS-Kg 
MODERATOR  MASS-Kg 
STRUCTURE  MASS-Kg 
REFLECTOR  MASS-Kg 
PRESSURE  VESSEL  MASS-Kg 
MISCELLANEOUS  MASS-Kg 
TOTAL  REACTOR  MASS-Kg 
NEUTRON  SHIELD  THICKNESS-m 
NEUTRON  SHIELD  MASS-Kg 
GAMMA  SHIELD  THICKNESS-m 
GAMMA  SHIELD  MASS-Kg 
TOTAL  SHIELD  MASS-Kg 


RADIATOR  PARAMETERS 
TOTAL  AREA-m2 
TOTAL  WEIGHT-Kg 


POWER  CONVERSION  PARAMETERS 
TURBINE  WEIGHT-Kg 
TURBINE  SPEED-RPM 
TURBINE  EFFICIENCY 
HEAT  EXCHANGER  WEIGHT-kg 
VAPOR  SEPARATOR  WEIGHT-kg 
CONDENSER  WEIGHT-kg 
GENERATOR  WEIGHT-Kg 
POWER  CONDITIONING  WEIGHT-Kg 
GEN  &  PC  RADIATOR  WEIGHT-Kg 
GEN  &  PC  RADIATOR  AREA-m2 
TOTAL  WEIGHT-Kg 


WEIGHT  AND  COST  SUMMARY 
REACTOR  WEIGHT-Kg 
SHIELD  WEIGHT-Kg 
POWER  CONVERSION  WEIGHT-Kg 
RADIATOR  WEIGHT-Kg 
MISCELLANEOUS  WEIGHT-Kg 
TOTAL  WEIGHT-Kg 


REACTOR+SHIELD  COST-M$ 
POWER  CONVERSION  COST-M$ 
RADIATOR  COST-M$ 

LAUNCH  COST-M$ 

TOTAL  COST-M$ 

Execution  terminated  :  0 


=  140.95736694 
=  165.93653870 
=  175.27770996 
=  316.23507690 
=  150000.00000000 
=  298.97799683 
=  2434.53515625 
=  2434.53515625 
=  O.OOOOOOOOE-01 
=  2374.32177734 
=  1392.97717285 
=  256.64099121 
=  2434.53515625 
=  8893.01074219 
=  O.OOOOOOOOE-01 
=  O.OOOOOOOOE-01 
=  6.23262003E-02 
=  4829.31933594 
=  4829.31933594 


2759.63330078 

17247.70703125 


5798.25781250 
10000.00000000 
0,. 78950000 
O.OOOOOOOOE-01 
1263.01953125 
1263.01953125 
1052.63159180 
2000.00000000 
2297.35815430 
367.57730103 
13674.28613281 


8893.01074219 

4829.31933594 

13674.28613281 

17247.70703125 

4464.43261719 

49108.75781250 


60.37825012 

40.26596069 

3.44954133 

62.85920715 

166.95295715 


TABLE  12 
Ranklne  System 

Mass  and  Cost  Algorithm  Summary 


Component _ _ Weight _ _ Cost 


Reactor  and  Shield 

-p - _ 

■  Marshall's  algorithm 
(Marshall,  1986) 

$4400/kg 

(estimate) 

Turbine 

See  Appendix  C 

$ 2000/kg 
(Gerry,  1985) 

Alternator 

0.1  kg/kW  (Gerry,  1985) 

$ 3000/kg 
(Gerry,  1985) 

Power  Conditioning 

0.2  kg/kW  (estimate) 

$10, 000/kg 
(estimate) 

Radiator 

- j 

12  kg/m2  temp  >  1000  K 

8  kg/m2  temp  650  to  1000  K 

5  kg/m2  temp  <  650  K 
(NASA  estimates) 

multiply  by  1.25  for  meteoroid 
losses  (estimate) 

See  Appendix  F  for  area 
calculation 

$ 200/kg 
(estimate) 

PC  6c  Alternator 

radiator 

5  kg/m^ 

(NASA  estimate) 

j  add  25%  for  meteoroid  losses 

$ 200/kg 
(estimate) 

Heat  Exchanger 

Preheater  6c  Boiler: 

.53[(hP  +  hb)m/1000] -74 

Superheater: 

14(hsm/1000) -7 
(NASA  estimates) 

hp  =  preheat  enthalpy  (J/kg) 
hb  ==  boiler  enthalpy  (J/kg) 
hg  «  superheat  enthalpy  (J/kg) 

m  =  flow  rate  (kg/s) 

$2000/kg 
(estimate- -same 
as  turbine) 

-51- 


TABLE  12  (cont.) 

Rankine  System 

Mass  and  Cost  Algorithm  Summary 
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THERMIONIC  CONTINUOUS  POffER  SYSTEM 


Like  th^  Brayton  and  Rankine  systems,  the  thermionic  system  provides 
continuous  power  and  is  powered  by  a  nuclear  reactor.  We  assume  the  use 
of  an  in-core  thermionic  reactor;  that  is,  the  thermionic  converters 
located  inside  the  reactor  are  cylindrical  in  geometry  and  enclose 
nuclear  fuel.  A  thermionic  reactor  will  contain  many  thermionic 
converter  elements .  Nuclear  fuel  heats  emitter  surfaces  which  causes 
electrons  to  "jump"  a  gap  between  the  emitter  and  collector.  The 
collector  collects  the  electrons  which  then  pass  through  an  electrical 
circuit  to  a  load  then  back  to  the  emitter.  The  collector  surfaces  are 
cooled  by  a  mixture  of  sodium  and  potassium  (NaK)  which  is  in  turn 
cooled  by  a  space  radiator.  In  terms  of  the  number  of  major  components, 
this  is  the  simplest  system.  It  consists  of  a  reactor,  power 
conditioning,  and  a  radiator. 


Our  model  for  the  performance  of  this  system  depends  heavily  on  Angrist 
(1982).  The  saturation  current  from  a  thermionic  emitter  can  be 
estimated  using  equation  52: 


AiT 


iJ-e 


2^  -^eATe 


and  the  back  current  from  the  collector  can  be  estimated  using 
equation  53: 


(52) 


where 


3,  -  We 


(53) 


jo  is  the  emitter  saturation  current  (A/m^) , 
jo  is  the  back  saturation  current  (A/m^)  , 

Ai  is  a  constant  (1,2  x  10®  a/m^) , 

To”  is  emitter  temperature  (K) , 

Tg  is  collector  temperature  (K) , 

</>g  is  the  emitter's  work  function  (eV)  , 

(i>c  is  the  collector's  work  function  (eV) , 

k  is  Boltzmann's  constant  (1.38  x  lO"^®  J/K,  but  it  must  be 
divided  by  1.6  x  IQ-i®  J/eV  to  get  8 . 62  x  lO'®  eV/K) ,  and 
AV  is  the  gap  voltage  drop  (eV) . 


Angrist  shows  that  the  minimum  voltage  drop  is  0.45  eV  and  that 
associated  with  it  is  a  saturation  current  ratio  of  0.4.  That  is, 
current  is  0.4  times  saturation  current. 


With  this,  the  power  per  unit  area,  P,  from  a  device  is: 

p  =  -^(je  -  jc)(<ie  -  4  -  AV)  .  (54) 

We  have  assumed  that  <f>e  -  <l>c  ~  AV  is  the  useful  voltage  thus  neglecting 
circuit  and  other  voltage  drops. 

To  calculate  the  efficiency  of  the  device,  we  need  to  do  an  energy 
balance  on  the  emitter.  Energy  carried  away  by  electrons,  Qg  (also 
called  electron  cooling),  is  given  as  follows: 
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Qe  -  je(^e  +  2kTJ  - 


(55) 


Radiation  losses,  Qr,  are; 

1 

- a  (Te-^  -  T,*) 

1  1 
—  + - 1 


(56) 


£q  is  emitter  thermal  emittance  (Angrist  suggests  0.3), 
Cg  is  collector  emittance  (Angrist  suggests  0.1),  and 
a  is  the  Stef  an- Boltzmann  constant  (5.7  x  10"®  W/m^k'')  . 


(57) 


Conduction  losses,  Q^,  are: 
k 

Qc  •=  -  (Te  -  Tc) 

g 

where  k  is  the  thermal  conductivity  of  cessium  in  the  gap  (Angrist 
suggests  .0124  W/mk) ,  and 
g  is  the  gap  thickness  (m) . 

Thus ,  the  heat  supplied  to  the  emitter  must  be  Q.  +  Qr  +  Qc 
efficiency  of  the  device,  tj,  is  as  follows: 


V  -  P/(Qe  +  Qr  +  Qc) 


(58) 


The  thermal  power,  P^hi  that  must  be  generated  by  the  reactor  is  given 
by  equation  59 : 


Pth  -  Pe/(»?»?pc) 

where  P.  is  the  electrical  power  required  by  the  load,  and 
ripe  is  power  conditioning  efficiency  (we  use  .95). 


(59) 


The  heat  which  must  be  dissipated  by  the  radiator  is  Qrad^ 


Qrad  -  Pthd  -  »?)  • 

We  assume  that  the  coolant  flow,  rate  through  the  reactor,  m,  is  5  kg/s 
per  MWth  of  heat  generated,  and  that  the  collector  temperature  is  equal 
to  the  reactor  coolant  outlet  temperature ,  Tj-^ .  The  reactor  inlet 
temperature,  Tj-i,  can  be  calculated  using  equation  61: 


Tri  ”  T^o  -  Pth/mCp 

where  Cp  is  the  NaK  coolant's  specific  heat  (950  J/kgK) . 


(61) 
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These  algorithms  were  written  into  a  computer  program  (TICYC)  using 
FORTRAN/77.  Collector  temperature  is  optimized  to  get  a  minimum  mass 
system.  As  collector  temperature  decreases,  efficiency  increases  and 
the  quantity  of  heat  that  must  be  rejected  decreases,  but  radiator 
temperature  decreases  and  the  radiator  rejects  less  heat  per  unit  area. 
Thus,  there  is  an  optimtam  collector  temperature. 

Table  13  lists  a  tjrpical  input  sequence  for  program  TICYC.  The  user 
must  specify  a  power  level,  an  operation  time,  and  then  whether  the 
program  is  to  minimize  mass  or  cost.  Following  tj^at,  default  parameter 
values  are  listed  and  the  user  may  change  them.  Table  14  shows  a 
reactor  inlet  temperature  (collector  temperature)  optimization  first, 
then  a  listing  of  parameters  for  the  optimized  system. 

The  program  estimates  component  mass  values.  Reactor  and  shield  masses 
are  estimated  using  algorithms  constructed  by  Marshall  (Ref.  1),  The 
mass  includes  fuel,  moderator  (if  any),  structure,  reflector,  gamma 
shield,  neutron  shield,  and  miscellaneous  mass.  Power  conditioning  is 
assumed  to  weigh  0.2  kg/kw  but  this  value  should  depend  strongly  on  load 
characteristics.  The  power  conditioning  is  assumed  to  have  a  95% 
efficiency  and  is  cooled  using  a  500  K  isothermal  radiator.  The 
system's  main  radiator  mass  is  summarized  in  Table  15  and  its  area  is 
calculated  using  the  algorithms  in  Appendix  F.  A  miscellaneous  mass  of 
10%  is  added  for  piping  and  structure. 

Costs  are  also  estimated,  but  they  are  not  accurate  enough  to  be  taken 
seriously  at  this  time. 
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TABLE  13 


THIS  PROGRAM  MODELS  A  THERMIONIC  SPACE  POWER 
SYSTEM.  THE  SYSTEM  CONSISTS  OF  A  LIQUID  METAL 
COOLED  THERMIONIC  REACTOR, 

A  RADIATOR,  AND  POWER  CONDITIONING. 


ENTER  VALUES  FOR  ELECTRICAL  POWER  IN  MW  AND 
OPERATING  TIME  IN  HOURS.  10  87600 

ENTER  "WEIGHT”  OR  "COST"  AS  PARAMETER  TO  BE 
MINIMIZED  (IN  CAPITAL  LETTERS) .  WEIGHT 

THE  FOLLOWING  ARE  DEFAULT  PARAMETER  VALUES. 
THEY  MAY  BE  CHANGED  IF  YOU  WISH. 


CYCLE  PARAMETERS : 

1.  1800.00000000 

2.  2.90000010 

3.  1.50000000 

4.  5.00000024E-04 

5.  8. 10000002E-02 

6.  0.20000000 

7.  0.94999999 


REACTOR 

PARAMETERS ; 

8. 

0.93000001 

9. 

18.00000000 

10. 

0.27000001 

11. 

24000.00000000 

12. 

0.10000000 

13. 

O.OOOOOOOOE-01 

14. 

O.OOOOOOOOE-01 

15. 

1.00000000 

16. 

4.99999992E+16 

17. 

25.00000000 

18. 

15.00000000 

19. 

2.00000000 

20. 

l.OOOOOOOOE+07 

EMITTER  TEMPERATURE,  K 

EMITTER  WORK  FUNCTION,  eV 

COLLECTOR  WORK  FUNCTION,  eV 

GAP  DIMENSION,  m 

GAP  EMITTANCE 

REACTOR  PRESSURE,  MPa 

POWER  CONDITIONING  EFFICIENCY 


FRACTIONAL  FUEL  ENRICHMENT 
CRITICAL  COMPACT  MASS,  Kg 
FUEL  +  MODERATOR  VOL  FRACT 
FUEL  POWER  DENSITY,  W/kg 
FUEL  BURNUP  FRACTION  LIMIT 
MODERATOR-TO-FUEL  RATIO 
MODERATOR  MOLECULAR  WEIGHT 
MODERATOR  DENSITY,  Kg/m3 
ALLOWED  PAYLOAD  NEUTRON  DOSE,  nvt 
PAYLOAD  SEPARATION  DISTANCE,  m 
PROTECTION  CONE  HALF  ANGLE,  DEG 
NEUTRON  SHIELD  MATL-B4C=1,  LIH=2 
ALLOWED  PAYLOAD  GAMMA  DOSE,  R 


RADIATOR 

21. 


PARAMETERS ; 

0.88000000 


RADIATOR  EMITTANCE 


ENTER  THE  NUMBER  OF  PARAMETERS  YOU  WISH  TO  CHANGE.  _0 
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TABLE  14 


OUTLET 

INLET 

TI 

RCT+SHLD 

POW  CONV 

RADIATOR 

TOTAL 

TOTAL 

TEMP  (K) 

TEMP  (K) 

EFFIC 

WT  (Kg) 

WT  (Kg) 

WT  (Kg) 

WT  (Kg) 

COST  (M$) 

1250.000 

1039.474 

0.026 

136720. 

3119. 

101993. 

266015. 

983. 

1200.000 

989.474 

0.091 

43689. 

3119. 

31211. 

85822. 

329. 

1150.000 

939.474 

0.108 

37355. 

3119. 

28794. 

76195. 

288. 

1100.000 

889.474 

0.113 

36021. 

3119. 

31035. 

77193. 

284. 

1150.000 

939.474 

0.108 

37355. 

3119. 

28794. 

76195. 

288. 

CYCLE  PARAMETERS 
EFFICIENCY 
THERMAL  POWER-MW 
MASS  FLOW  RATE-Kg/s 
REACTOR  OUTLET  TEMP-K 
REACTOR  INLET  TEMP-K 
REACTOR  PRESSURE-MPa 


0.10808214 

97.39181519 

486.95907593 

1150.00000000 

939.47369385 

0.20000000 


REACTOR  PARAMETERS 
BURNUP  MASS-Kg 
INITIAL  CRITICAL  MASS~Kg 
END  CRITICAL  MASS-Kg 
TOTAL  BRNUP+CRIT  MASS-Kg 
SPECIFIC  POWER-W/Kg 
MASS  FOR  SPECIFIC  POWER  LIM-Kg 
MASS  FOR  ALLOWED  BURNUP-Kg 
FUEL  MASS-Kg 
MODERATOR  MASS-Kg 
STRUCTURE  MASS-Kg 
REFLECTOR  MASS-Kg 
PRESSURE  VESSEL  MASS-Kg 
MISCELLANEOUS  MASS-Kg 
TOTAL  REACTOR  MASS-Kg 
NEUTRON  SHIELD  THICKNESS-m 
NEUTRON  SHIELD  MASS-Kg 
GAMMA  SHIELD  THICKNESS-m 
GAMMA  SHIELD  MASS-Kg 
TOTAL  SHIELD  MASS-Kg 


397.94503784 

235.27148437 

260.43603516 

658.38110352 

24000.00000000 

5275.38964844 

4811.15527344 

5275.38964844 

O.OOOOOOOOE-01 

11955.35156250 

2722.29345703 

309.63922119 

5275.38964844 

25538.06445312 

0.19794284 

1241.26269531 

7.31635392E-02 

10575.74707031 

11817.00976562 


RADIATOR  PARAMETERS 

RADIATOR  INLET  TEMPERATURE,  K 
RADIATOR  OUTLET  TEMPERATURE,  K 
HIGH  TEMPERATURE  AREA,  m2 
MEDIUM  TEMPERATURE  AREA,  m2 
LOW  TEMPERATURE  AREA,  m2 
TOTAL  AREA-m2 
TOTAL  WEIGHT-kg 


1110.00000000 

899.47369385 

829.06091309 

1155.91284180 

O.OOOOOOOOE-01 

1984.97375488 

28794.05078125 


POWER  CONVERSION  PARAMETERS 

POWER  CONDITIONING  WEIGHT-kg 
PC  RADIATOR  WEIGHT-kg 
TOTAL  WEIGHT-Kg 


2000.00000000 

1119.22570801 

3119.22583008 
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TABLE  14  (cont.) 


WEIGHT  AND  COST  SUMMARY 
REACTOR  WEIGHT-Kg 
SHIELD  WEIGHT-Kg 
POWER  CONVERSION  WEIGHT-Kg 
RADIATOR  WEIGHT-Kg 
MISCELLANEOUS  WEIGHT-Kg 
TOTAL  WEIGHT-Kg 


25538.06445312 

11817.00976562 

3119.22583008 

28794.05078125 

6926.83496094 

76195.18750000 


REACTOR+SHIELD  COST-M$ 
POWER  CONVERSION  COST-M$ 
RADIATOR  COST-M$ 

LAUNCH  COST-M$ 

TOTAL  COST-M$ 

Execution  terminated  :  0 


164.36231995 

20.22384453 

5.75881004 

97.52983856 

287.87481689 
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TABLE  15 

Thermionic  System 
Mass  and  Cost  Algorithm  Summary 


Component _ _ _ _ Weight _ _ Cost 


Reactor 

Marshall's  algorithm 
(Marshall,  1986) 

$4400Ag 

(estimate) 

Power  Conditioning 

0.2  kg/kW  (estimate) 

$10, 000/kg 
(estimate) 

PC  Radiator 

5  kg/m^ 

$200Ag 

(NASA  estimate) 

multiply  by  1.25  for  meteoroid 
losses  (estimate) 

(estimate) 

Radiator 

12  kg/m^  temp  >1000  K 

$200/kg 

8  kg/m^  temp  650  to  1000  K 

5  kg/m^  temp  <650  K 
(NASA  estimates) 

multiply  by  1.25  for  meteoroid 
losses  (estimate) 

muliply  by  1.2  for  heat  exchanger 
(estimate) 

(estimate) 
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SUMMARY  AND  CONCLUSIONS 


This  docximent  has  described  five  multimegawatt  space  power  system  models: 
Gas  cooled  reactor  powered  Brayton  cycle, 

Liquid  metal  cooled  reactor  powered  Rankine  cycle, 

Liquid  metal  cooled  reactor  thermionic  system. 

Gas  fooled  reactor  powered  open  turbine  generator,  and 
Hydrogen- oxygen  combustion  open  turbine  generator. 

The  mathematical  algorithms  used  to  estimate  system  performance  and  mass 
have  been  described.  These  models  have  formed  the  foundation  of  our 
efforts  to  evaluate  the  many  proposed  multimegawatt  space  power  systems . 

Showing  some  examples  of  how  these  models  have  been  used  would  be  helpful 
here,  but  as  an  alternative,  we  refer  the  reader  to  other  doctiments  which 
describe  system  studies  conducted  using  the  models  (Edenburn  1988-1  and 
1988-2,  and  Edenburn  1990). 
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APPENDIX  A 


Specific  Power  Limit  For  A 
Gas  Cooled  Particle  Bed  Reactor 

The  specific  power  extracted  from  a  reactor  core  can  be  limited  by 
either  fuel  temperature  or  pressure  drop  considerations.  The  specific 
power  limit  determined  by  temperature  is  the  maximum  specific  thermal 
power  (power  divided  by  fuel  mass)  that  can  be  removed  without  causing 
the  fuel's  temperature  limit  to  be  exceeded.  To  determine  this  specific 
power  limit  we  will  consider  a  fuel  particle  comprised  of  a  central  fuel 
kernel  which  is  coated  by  layers  of  other  materials. 


The  temperature  at  the  center  of  the  fuel  kernel  can  be  found  using  a 
standard  radial  heat  conduction  analysis. 


1  q  rl^  1  q  1  n  1 

Tf  -  Tc  + - + - qri^  s  — 

3  h  rn2  6  ki  3  i=2  ki 


1 


1 


ri-1  ri  J 


(A-1) 


where  Tf  is  the  temperature  at  the  center  of  the  fuel  kernel  (K) , 

Tq  is  coolant  temperature  (K) , 

q  is  the  volumetric  heat  generation  in  the  fuel  kernel  (W/m^) , 
h  is  the  coolant's  convection  coefficient  (W/m^  K) , 
rj_  is  the  radius  of  the  i^^  layer  (m)  (the  fuel  kernel  is 
number  1),  and 

ki  is  the  thermal  conductivity  of  the  i^^^  layer  (W/m  K)  . 


The  term  following  T^  is  the  temperature  rise  from  the  coolant  to  the 
surface  of  the  particle.  Following  that  is  the  term  representing  the 
temperature  rise  from  the  edge  to  the  center  of  the  fuel  kernel.  The 
summation  term  is  the  temperature  rise  for  all  of  the  intermediate 
layers.  If  we  let  5  represent  the  kernel  and  layer  temperature  rise 
terms  (divided  by  q) ,  let  R  be  the  ratio  of  kernel  to  particle  radius 
and  D  be  the  particle  diameter  Equation  (A-1)  can  be  abbreviated. 


Tf  =  Tc  +  I  J  r3d  +  5q 


(A-2) 


A-1 


Solving  this  for  q,  we  get  equation  A- 3. 


(A- 3) 


Since  q  is  the  volumetric  heat  generation,  q/^f  is  the  fuel's  specific 
power,  is  the  density  of  the  fuel  kernel.  Its  limiting  value  occurs 

when  Tg  is  at  its  maximum  -  the  reactor's  outlet  temperature,  Tout  " 
when  the  fuel  temperature  is  at  its  maximum  value,  Tfn,ax- 


Ps 


*Ima3C 


'•fmax 


6h 


'Pput _ 

+  S^t 


(A-4) 


Pg  is  the  specific  power  limit  due  to  temperature  considerations. 
However,  we  have  not  yet  evaluated  h.  Before  evaluating  h,  we  will 
define  the  rest  of  the  terms  to  be  used  in  the  thermal  analysis: 


A  is  the  bed's  cross  sectional  area  (m^)  , 

Ap  is  the  total  particle  surface  area  (m^) , 

Cp  is  the  coolant's  specific  heat  (J/kgK), 

D  is  the  fuel  particle  diameter  (m) , 

h  is  the  heat  transfer  coefficient  between  the  fuel  particles  and 
the  coolant  (W/m^K)  , 

k  is  the  coolant's  thermal  conductivity  (W/mK) , 

L  is  the  length  of  fuel  bed  through  which  the  coolant  flows  (m) , 

m  is  the  coolant  flow  rate  (kg/s), 

Pr  is  the  coolant's  Prandtl  number, 

R  is  the  ratio  of  fuel  kernel  radius  to  fuel  particle  radius , 

Re  is  the  coolant's  Reynolds  number. 

Tin  is  the  coolant  inlet  temperature  (K) , 

T^„  is  the  maximum  allowed  fuel  particle  surface  temperature  (K) , 
Tout  is  the  coolant  outlet  temperature  (K) , 

U  is  the  coolant's  free  stream  velocity  (m/s), 

Vj,  is  the  total  fuel  bed  volxime  (m®)  , 

Vf  is  the  total  fuel  kernel  volume  (m^)  (a  fuel  particle  consists 

of  a  UC  fuel  kernel  surrounded  by  one  or  more  carbon  coatings) , 
Vp  is  the  total  fuel  particle  volume  (m^)  , 

^  is  the  coolant  density  (kg/m®)  , 

is  the  fuel  kernel  density  (kg/m®)  ,  and 
fi  is  the  coolant's  viscosity  (kg/ms) . 

Eckert  (1959),  has  an  expression  for  h  for  a  packed  bed  when  Re  >  500. 


h 


J  Up£  jo.7  333 


(A-5) 


The  term  in  parenthesis  is  a  Reynolds  number  and  is  generally  greater 
than  500  for  our  applications. 


A-2 


We  need  to  find  U,  the  coolant's  free  stream  velocity  to  evaluate  h  so 
we  can  determine  the  value  of  Pg.  To  find  U,  we  will  use  the  energy 
balance  equation  for  the  reactor,  Equation  (A-6). 

p  „  mCpCTput  -  Tin)  (A-6) 

s  PfVf 

We  can  put  this  in  a  useful  form  using  the  following  relations: 


m  =  PHA,  (A-7) 
A  -  Vb/L,  (A- 8) 
Vf  =  R^Vp,  and  (A- 9) 
Vp  =  .65Vb.  (A- 10) 


This  last  equation  assumes  that  the  particle  packing  density  is  65%. 
(A-6)  can  now  be  written  as, 

U^Cp(Tout  -  Tin) 

Ps - -  .  (A-11) 

.eSPfR^L 

We  can  solve  this  for  U,  use  this  value  of  U  to  find  h,  and  substitute 
the  value  of  h  into  Equation  (A-4)  to  evaluate  Pg. 

Letting 


ai  = 


.282 


k  P 


0.33 


Cp/i(Tout  “  "^in) 


0.7 


(A-12) 


Equation  (A-4)  becomes: 


T  -  T 
fmax  out 


■0:7 -^'"f 


(A-13) 


or 


a,P  +  SP f?s  =  T.  -  T  ^  .  (A-14) 

Is  ^  ^  fmax  out 

This  equation  can  be  solved  for  Pg. 

That  gives  us  the  specific  power  l.imit  associated  with  fuel  temperature. 
Now  we  will  turn  to  the  specific  power  limit  imposed  by  pressure  drop 
considerations.  The  pressure  drop  across  the  particle  bed  is  given  by 
Equation  (A-15)  (Eckert  and  Drake,  1959). 


A-3 


+  1.75 


(A-15) 
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where  c  is  the  void  fraction  of  the  bed  (1-0.65  «=  0.35);  see  Equation 
(A-10), 

P  is  the  outlet  pressure  from  the  reactor  (Pa) > 

A  is  the  allowed  fractional  pressure  drop,  and 
bars  over  parameters  denote  mean  values. 

U  will  vary  through  the  bed  because  density,  P,  changes  with 
temperature,  however  is  constant  (this  is  not  precisely  true  for 
radial  flow  but  is  a  reasonable  approximation  if  the  inside  and  outside 
fuel  bed  diameters  are  about  the  same) .  The  coolant  can  be  treated  as 
an  ideal  gas,  thus  we  can  evaluate  U  as  in  Equation  (A-16) . 


u-f 


where  U  is  the  coolant  velocity  (m/s)  at  the  bed's  exit. 


(A-16) 


We  have  also  estimated  the  pressure  drop  through  flow  channels  from 
Powell  (1985). 


AP  =  2.9  X  10^  Tin/Tout  +  7.6  x  10^  j  (moderated) 

AP  =  6  X  10^  TinAout  +  7.6  x  10^  j  (unmoderated) 

The  two  pressure  drops  are  added  to  get  the  total  pressure  drop. 


(A-17) 


AP  =  740  1  +  ^+13.3  1  +  +  rPU^, 

out  D  out  ^ 


(A-18) 


where  F  is  defined  by  the  bracketed  parts  in  Equation  (A-17). 


Again  the  energy  balance  Equation  (A-11)  must  be  satisfied. 

Substituting  the  value  of  U  from  the  energy  balance  equation  into  (A-18) 
gives  us  an  expression  for  the  specific  power  limit  based  on  the  allowed 
pressure  drop. 


AP  -  a.P  +  a.P 
2  s  3  s 


(A-19) 


a2  =  481  I  1  + 


^in  1/*^'  V' 
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For  specified  values  of  A  and  P  we  can  find  the  value  of  Ps. 

The  two  values  of  Pg  will  not  be  the  same  and  the  smaller  of  the  two 
should  be  selected. 


A  summary  of  coolant  properties  is  given  in  Table  A-1. 


A- 5 
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TABLE  A-1 


Coolant  Properties 


ProDertv 

H9 

He 

>(i) 

f  T  ^  nO.43 

0  49  1  1 

•  [  1200  J 

f  T  ^  >0.4 
0.30  [-li] 

■(S) 

23x10-6  1  ] 

^  [  1200  J 

r  T  ^  >0.4 
^1x10  ^  ^22] 

15,400 

5188 

P 

P 

\^out 

R^T  ^ 

G  out 

Pr 

0.72 

0.72 

"c  (  EgK  , 

4126 

2077 

V 

1.41 
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APPENDIX  B 


WRAPAROUND  REACTOR  SHIELD  SIZE  CALCULATION 

Components  behind  the  reactor  are  protected  from  radiation  by  a  shadow 
shield.  In  Marshall's  model  (Marshall,  1986),  the  shield  is  shaped  like 
a  disc  whose  radius  depends  on  the  cone  half- angle  of  protection.  For 
cone  half -angles  above  a  certain  value,  a  "wraparound"  shield  that 
protects  a  cone  half- angle  of  90°  would  be  lighter  than  a  disc  shield. 


REACTOR 


Marshall's  model  calculates  the  required  thickness  of  alternating 
neutron  and  gamma  shield  layers  ti,  t2,  t3,  and  ti^.  The  wraparound 
shield  calculation  uses  these  thicknesses  to  calculate  the  weight  of  a 
disc  bottom  shield  and  a  tapered  side  or  "wraparound"  shield  as  shown  in 
Figure  B-1.  Lg  is  the  reactor's  length  and  r+T  is  its  radius.  The 
weights  of  the  four  layers  are  given  as  Wi,  W2,  W3,  and  W4. 

PyQ  and  P^iQ  are  the  densities  of  the  gamma  and  neutron  shields 
respecively , 

If  the  weight  of  the  wraparound  shield  is  less  than  that  of  the  disc 
shadow  shield,  its  value  is  used  as  the  shield  weight. 
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APPENDIX  C 
Turbine  Models 


The  gas  turbine  model  described  in  this  appendix  was  formulated  by 
Hudson  (1988).  Modifications  to  his  original  formulation  adjust 
material  properties  in  each  stage  for  temperature  and  add  algorithms  for 
blade  cooling. 


An  axial  flow  gas  turbine  is  composed  of  one  or  more  stages  which 
extract  energy  from  a  working  fluid.  Fixed  nozzles  accelerate  the 
working  fluid  and  direct  it  onto  rotating  blades  which  produce 
mechanical  work.  We  determine  the  turbine's  speed,  disk  (to  which  the 
blades  are  attached)  radius ,  and  blade  length  with  help  from  four 
equations:  1)  the  blade  stress  equation  (equation  C-1);  2)  the  disk 
stress  equation  (equation  C-2);  3)  the  working  fluid  continuity  equation 
(equation  C-3);  and  4)  the  energy  equation  (equation  C-4) . 


Ob  =  (R  +  L/2)w^ 
f 


1  „2  2 

-79  =  R  w 


A  “  pV^27r  (R  +  ^)L  ,  and 
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1  fv?  -  1  +  C  (T.  -  T  ) 

2  ^  in  outj  p  in  out' 


Power 


We  will  also  need  the  isentropic  expansion  relation, 
state,  and  an  efficiency  algorithm. 


(C-1) 

(C-2) 

(C-3) 

(C-4) 

an  equation-of- 


where 

CTb  is  blade  specific  strength  (strength  divided  by  density), 
Tf  is  a  taper  factor  (.7), 

R  is  disk  radius, 

L  is  blade  length, 
w  is  turbine  speed  in  rad/s, 

Ojj  is  disk  specific  strength, 

m  is  mass  flow  rate, 
p  is  the  working  fluid's  density, 

Va  is  the  working  fluid's  axial  velocity, 

Vin  is  the  working  fluid's  speed  entering  a  stage, 

Vout  is  the  working  fluid's  speed  leaving  the  stage. 

Tin  Tout  are  corresponding  temperatures ,  and 
Cp  is  the  fluid's  specific  heat. 


Equation  C-1  sets  the  blade's  root  tensile  stress  at  its  maximum  allowed 
strength  and  equation  C-2  sets  the  disk's  greatest  tensile  stress  at  its 
maximum  allowed  strength. 
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Equation  C-3  requires  that  the  mass  flow  rate  through  a  turbine  stage  be 
equal  to  the  fluid's  density  multipled  by  its  volume  flow  rate.  V^,  in 
this  equation,  can  be  written  in  terms  of  the  turbine's  blade  speed,  U, 
by  defining  the  turbine's  work  coefficient,  4>. 

V  +  V 

Tin  Tout  (C-5) 


where 


U-(R+'l/2)«  , 

and  by  specifying  that  it  is  to  be  a  50%  reaction  turbine  with  constant 
axial  velocity  as  in  equation  C-7. 

V  V  «»  U  (C"7) 

Viin  Viout  ^ 

This  specification  results  in  half  of  the  stage's  enthalpy  drop  occuring 
in  the  nozzles  and  half  in  the  blades.  A  100%  reaction  stage  would 

extract  all  of  the  enthalpy  drop  in  the  blades  and  an  impulse  stage  with 

0%  reaction  would  extract  all  in  the  nozzles . 

Viin  is  the  tangential  fluid  velocity  and  is  the  leaving  tangential 

velocity. 

Equations  C-5  and  C-7  can  be  combined  to  evaluate  Vxouf 

Vut  -  -  1)  ■ 

We  must  also  introduce  the  angle,  a,  at  which  the  fluid  leaves  the 
relative  to  the  rotor  blades  as  in  Figure  C-1. 


(C-8) 

rotor 


Vaout 


Figure  C-1.  Turbine  Velocity  Vector  Diagram 


Notice  that  U  was  added  to  Vj^^t  to  get  the  tangential  velocity  relative 
to  a  blade.  We  can  now  get  an  expression  for  V^  to  use  in  equation  C-3. 
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V 


aout 


=  (V, 


Tout 


+  U)/tan  a 


(C-9) 


V  =  U  t  ^  ^ 
aout  2  tan  a 


(C-10) 


Using  this,  the  ideal  gas  relation,  and  Equation  (C-6),  equation  (C-3), 
the  continuity  equation,  can  be  put  in  its  final  form. 


where 


riiR  T  ^tan  a 
S  out _ 

jr  P  (<f>  +  1) 
out 


Rg  is  the  fluid's  gas  constant, 

Tout  is  its  outlet  temperature  from  a  stage,  and 
^out  is  i^s  outlet  pressure. 


(C^ll) 


We  can  find  a  stage's  blade  length,  L,  disk  radius,  R,  and  angular 
speed  w,  by  simultaneously  solving  Equations  (C-1),  (C-2),  and  (C-11); 
however,  there  is  a  special  process  we  must  use.  First,  we  need  to  know 
the  turbine's  outlet  temperature  and  pressure.  To  find  these  we  use  a 
specified  turbine  pressure  ratio,  an  estimate  of  turbine  efficiency 
(Equation  C-12) ,  and  the  isentropic  expansion  relation,  modified  for  a 
turbine  efficiency  of  (Equation  C-13).  ryg  is  the  turbine's  stage 
efficiency  and  is  calculated  using  Equations  (C-22)  and  (C-23)  given 
later  in  this  appendix.  (These  equations  are  for  turbines  with  many 
stages . ) 

1  ^  j  ’?sRg/Cp 

T  ^  =  T.  {  I  -  ] 

out  in  I  t  't  p  J 

where  Tout  is  the  turbine's  outlet  temperature, 

Tin  IS  Its  inlet  temperature, 
is  the  estimated  efficiency, 

Rp  is  the  specified  pressure  ratio,  and 
Cp  is  the  working  fluid's  specific  heat. 

(Cp  will  be  constant  for  helitun,  but  for  other  gasses,  it  will 
vary  with  temperature  and  can  be  evaluated  using  the  enthalpy 
relations  in  Appendix  G.  When  Cp  is  not  constant,  an  iterative 
process  must  be  used  to  solve  Equation  (C-12),) 

The  turbine's  last  stage  will  be  its  largest  stage,  and  since  we  assume 
that  all  of  the  stages  must  have  the  same  angular  speed,  the  blades  in 
the  last  stage  will  experience  the  greatest  stress.  Because  of  this,  we 
start  our  process  with  the  last  stage,  using  the  turbine's  outlet 
temperature  and  pressure,  to  find  the  turbine's  speed.  Thus,  we 
calculate  w,  R,  and  L  for  the  last  stage.  These,  of  course,  depend  on 


(C-12) 

(C-13) 
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blade  and  disk  strength  which  are  functions  of  temperature.  Blade  and 
disk  temperatures  for  the  last  stage  are  presumed  to  be  equal  to  the 
working  fluid's  outlet  temperature  or  to  a  prescribed  temperature  if  the 
blades  are  cooled.  Material  strength  and  its  dependence  on  temperature 
are  catalogued  at  the  end  of  this  appendix.  The  value  of  w,  which  we 
have  just  calculated  for  the  last  stage,  is  the  value  of  w  we  will  use 
for  all  of  the  turbine's  stages  unless  we  have  prescribed  a  lower  value 

for  w. 

The  next  step  in  the  process  is  to  calculate  R  and  L  for  each  turbine 
stage  starting  with  the  first  stage.  We  calculate  R,  disk  radius,  using 
Equation  (C-2).  Disk  strength  properties  are  evaluated  at  the  stage's 
inlet  temperature  or  at  a  prescried  temperature  if  the  disk  and  blades 
are  cooled.  To  find  L  for  the  first  and  subsequent  stages,  we  must 
require  that  energy  and  working  fluid  be  conserved  as  working  fluid 
flows  through  a  stage.  Equation  (C-4)  is  the  energy  balance  equation. 
The  stage's  power  is  equal  to  the  force  exerted  on  the  blades  by  the 
fluid,  due  to  its  momentum  change,  multiplied  by  the  blade's  speed. 

Force  =  m(Vtin  +  ^^out)  •  (C-14) 


Power  =  m(Vtin  +  Vtout)^  •  '' 

Using  Equation  (C-5)  we  get 

Power  =  =  m<^(R  +  L/2)2w2  ,  (C-16) 

Vout  put  in  terms  of  U  and  finally  L  as  follows : 

V^=v^  +V^  =  [(^-l)U/2l^+  [(,^+l)U/(2tanQ)  1  .  (C-17) 

out  tout  aout  L  J  L  ^ 


These  came  from  Equations  (C-8)  and  (C-10). 

With  a  little  regrouping,  and  using  U  ^  (R  +  L/2)w,  the  energy  equation 
can  be  rewritten  as  in  (C- 18). 


^  <hn  -  ’'out*  - 


(C-18) 


r 

~  r\  O 

A  =  " 

^  ^ 

-  ir  ,  (<!>  + 1) 

2 

4  4tan  a 

J 

0) 


The  continuity  equation,  (C-11) ,  restated  here  is  as  follows: 


(C-19) 


”^g'’^out^^"  °  .  =  r  R  +  L  1  2 

.  P,ut(^  -^1)  I  2  J 


(C-11) 
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^in  Tin  known  from  the  outlet  conditions  of  the  previous  stage 
or,  for  the  first  stage,  from  the  turbine's  inlet  conditions,  R  and  w 
were  already  calculated  and  a  and  ^  are  specified  values ,  The  unknowns 
are  Tout*  T.,  and  Pout-  ^p  known  as  a  function  of  temperature.  To 
find  the  three  unknowns,  we  guess  values  for  Pout  ^p  (our  guess  is 
that  Pout  =  Pin  Cp  “  Cp(Tin))  solve  the  two  equations 
simultaneously  to  find  Tout  T**  using  an  iterative  procedure.  Then 
we  find  a  new  estimate  for  Pout  using  the  isentropic  expansion  relation 
and  the  stage's  efficiency,  and  we  estimate  a  new  Cp  based  on  the 
enthalpy  change  between  Ti^  and  Tout- 


P  ^  “  P. 
out  in 


in 

out 


Cp/R 


(C-20) 


ip  ^  _  iji 

out  in 


(T.  -  T  J 

s  in  out 


(C-21) 


f7s  is  stage  efficiency.  The  value  of  stage  efficiency  was  formulated  by 

•jif 

Steve  Hudson  (1987)  using  empirical  data  from  a  number  of  sources.  Tout 
is  the  stage's  outlet  temperature  for  an  ideal,  isentropic  expansion, 
and  Tout  is  the  real  outlet  temperature. 


ri  =  2.49C  -  1.623C^  .  (C-22) 

s 

C  -  2  sin  a/{<f>  +  1)  .  (C-23) 

The  new  values  of  Cp  and  Pout  used  in  the  energy  and  continuity 
equations  to  calculate  new  values  of  L  and  Tout*  iterative 

process  is  continued  until  it  converges. 

Notice  that  we  have  not  used  the  blade  strength  equation.  Since  blade 
stress  is  greatest  in  the  last  stage,  where  it  is  equal  to  the  blade's 
strength,  blade  stress  will  usually  be  less  than  strength  in  the  other 
stages.  But,  since  temperatures  change  from  stage  to  stage,  it  is 
possible  to  exceed  blade  strength  in  an  intermediate  stage.  We  check  to 
see  if  blade  strength  has  been  exceeded  using  Equation  (C-1) .  If  it 
has,  a  lower  value  for  w  must  be  specified. 

In  this  manner,  we  step  stage-by-stage  through  the  turbine  and  find 
blade  length  and  disk  radius.  We  stop  when  the  turbine's  outlet 
temperature  is  less  than  or  equal  to  the  desired  outlet  temperature  from 
Equation  (C-13) .  It  is  likely  that  we  have  overshot  the  desired  outlet 
temperature  substantially,  and  this  is  corrected  by  multiplying  the  disk 
radii  by  a  number  less  than  one  to  reduce  the  turbine's  power 
extraction.  We  iterate  the  disk  size  until  we  get  the  desired  outlet 
temperature . 
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Recall  that  we  found  the  desired  outlet  temperature  using  an  estimate 
for  turbine  efficiency  in  equation  C-12.  We  can  find  a  more  accurate 
turbine  efficiency  using  the  modified  isentropic  expansion  relation. 

“  ^"^in  ■  ^out  calculated)/ -  T^;^)  .  (C-24) 


T.  /(R  calculated)^^^P 
in  p 


(C-25) 


The  above  computations  give  us  the  disk  radius,  R,  and  blade  length,  L, 
for  each  stage.  We  find  the  mass  of  a  stage  by  assuming  that  the 
stage's  axial  length  depends  on  blade  length.  Aspect  ratio,  blade 
length  divided  by  sta^e  length,  is  given  by  equation  C-26,  and  stage 
length  is  given  by  C-27. 


Ap  =  4  -  3.2e"-^^ 


(C-26) 


Sl  -  L/Ar,  R/25,  or  .004  m  whichever  is  greater. 


(C-27) 


The  rotor  disk  fills  100%  of  its  associated  volume  plus  20%  of  the 
stator  volume  and  the  blades  fill  30%  of  their  associated  volume  for 
both  the  stator  and  rotor. 

disk  mass  -  120%7rR2p„SL  .  (C-28) 

blade  mass  -  30%jr[(R+L)2  .  r2]p^2Sl  .  (C-29) 

casing  mass  =  27r(R+L)t^n,2SL  •  (C-30) 

T}^  is  the  turbine  material's  density. 

Casing  thickness,  t,  is  found  using  a  hoop  strength  calculation. 


(R+L)P.  P  .  , 

t  -  - is  or  whichever  is  greater.  (C- 

m  b 

Recall  that  Brayton  cycle  efficiency  depends  on  blade  coolant  flow 
rate,  m.  A  stage's  blade  coolant  flow  rate,  liig,  depends  on  the 
following  parameters: 

h  -  the  convection  coefficient  between  the  working  fluid  and  the 
blade , 

A  -  total  blade  surface  area. 

Tin  -  stage  inlet  temperature, 

Tj,  -  maximum  allowed  blade  temperature , 

Cp  -  working  fluid  specific  heat,  and 

Tg  -  coolant  inlet  temperature,  (equal  to  compressor  outlet 
temperature  in  a  Brayton  cycle) . 


(C-31) 
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El  Wakil  (1984)  gives  a  graphic  relation  between  mgCp/hA  and 
what  he  calls  blade  cooling  effectiveness,  (Txn“'^b)/C^in"'^c)  » 
various  types  of  blade  cooling.  The  best  was  a  combination  of  film  and 
convective  cooling,  and,  for  this,  a  close  fit  to  his  data  is  given  by 
Equation  (C-32). 


m  C 
s  P 

hA 


=  4 


(C-32) 


Thus,  if  we  have  values  for  h  and  A,  we  can  calculate  a  value  for  ihg . 

El  Wakil  proposes  values  of  h  between  200  and  500  B/ft^hr  F  for  air.  In 
general,  for  laminar  flow  over  exterior  surfaces,  h  can  be  estimated 
using  Equations  (C-33)  and  (C-34)  since  Prandtl  number  is  fairly 
constant  among  gases. 

Nu  a  ,  and  (C-33) 


h  a  k 


Vp 

/iL 


(C-34) 


where  Nu  is  Nusselt  number. 

Re  is  Reynolds  number, 
k  is  gas  conductivity, 

V  is  gas  velocity 

p  is  gas  density, 

fi  is  gas  viscosity,  and 

L  is  the  dimension  of  the  surface. 


If  gas  velocity  and  blade  geometry  are  similar  for  turbines  using  two 
different  gases,  then  Equation  (C-34)  can  be  simplified  as  follows: 


h  a  k 


(C-35) 


Thus,  helium  should  have  1,4  times  the  convective  coefficient  of  air 

since  its  kj p//i  is  1.4  times  that  of  air  (k  is  thermal  conductivity, 
p  is  density,  and  /z  is  viscosity).  We  selected  400  B/ft2hr  F  for  air 
and  multiplied  by  1.4  to  get  an  h  of  560  B/hr  ft^  F  or  3000  J/m^sK  for 
helium.  For  mixtures  of  He  and  X^,  k  and  p  are  fairly  constant  until 
the  mass  fraction  of  Xe  gets  very  large,  greater  than  0.8.  Thus,  we 
assume  that  h  for  a  He-Xe  mixture  depends  on  density,  and,  since  density 
is  inversely  proportional  to  the  mixture's  gas  constant,  we  can  estimate 
h  for  the  mixture  as  follows : 

-  V  ■ 
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Thus, 


He-Xp 


3000 


J _ 

m^sk 


2079 

^e-Xe 


(C-37) 


Keeping  in  mind  that  each  blade  has  two  sides,  that  there  are  two  sets 
of  blades  (nozzles  and  rotor  blades)  in  each  stage,  and  assuming  that 
the  spacing  between  blades  is  (1/3)L,  we  can  find  blade  surface  area,  A. 

A  -  247rRSL  .  (C-38) 


From  Equations  (C-32),  (C-36),  and  (C-38),  we  calculate  mg  for  each 
stage  and 

add  each  of  them  to  find  the  total  coolant  flow  rate,  m.  When  the 
maximum  allowed  blade  temperature  is  above  the  stage  inlet  temperature, 
no  coolant  is  necessary. 
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Turbine  Material  Strength  Catalog 


We  assume  that  if  the  blades  are  being  cooled,  the  disk  is  maintained  at 
the  inlet  coolant  temperature. 

These  specific  strength  (strength  divided  by  density)  values  are  fits  to 
data  compiled  by  Steve  Hudson  and  represent  values  that  result  in  1% 
creep  over  7  years  except  for  carbon  composites  where  creep  data  is  not 
available.  For  carbon  composites,  the  values  are  estimates  based  on 


ultimate  tensile  strength. 

Ni  Superallov: 

below  875  K  a  =  110  kJ/kg 

875  to  1350  K  a  =  3.961  +  5 . 377xlO5e-0  00975T  kJ/kg  (C-39) 

above  1350  K  -  unusable 

=  8500  kg/m3  (C-40) 

TZM  Mo  Refractory: 

below  1000  K  a  =  40  kJ/kg 

1000  to  1800  K  a  -  -2.21  +  1.046xl03e~o.oo32iT  kJ/kg  (C-41) 

above  1800  K  -  unusable 

==  10,200  kg/m3  (C-42) 

W-HfC  W  Refractory: 

below  1375  K  a  =  36  kJ/kg 

1375  to  2000  k  a  -  -67.2  +  2 . 567xl02e“0  000663T  kJ/kg  (C-43) 

above  2000  K  -  unusable 

=  19,300  kg/m3  (C-44) 

SIqN/,  Ceramic: 

below  870  K  a  124  kJ/kg 

870  to  1590  k  a  =  212  -  .1014  T  (C-45) 

above  1590  K  -  unusable 

-  3,000  kg/jn3  (C-46) 

Carbon  Composite: 

all  temperatures  a  ^  167  kJ/kg  (C-^7) 

=  1800  kg/m3  (C-48) 
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I  have  used  the  preceding  turbine  model  to  develop  simple  turbine  mass 
algorithms  for  superalloy  and  carbon-carbon  composite  hydrogen- -oxygen 
combustion,  hydrogen  and  helium  turbines,  and  similar  algorithms  were 
developed  for  potassitim- vapor  turbines  by  Steve  Hudson. 

Equation  (C-49)  shows  the  general  form  of  the  simplified  mass 
algorithms,  however,  modifications  are  required  for  some  specific 
turbines.  Also,  we  at  times  used  different  sets  of  coefficients  for 
different  parameter  ranges. 

mass  (kg)  -  Mq  P«  Pi/^  Rp7  (1  +  uepTd) ,  (C-A9) 

If  Ti  <  Id  then  T^. 

The  parameters  are  described  as  follows: 

P  is  the  shaft  power  generated  in  MW, 
is  turbine  inlet  temperature  in  K, 
is  turbine  inlet  pressure  in  MPa, 

Rp  is  turbine  pressure  ratio, 

<l>  is  the  work  coefficient, 
u>  is  turbine  speed  in  RPM, 

Td  is  disk  temperature. 

Mo,  ctf  Pj  7,  and  p  are  coefficients  to  be  derived  for 

each  type  of  turbine. 

With  blade  and  disk  cooling,  Td  is  lower  than  Ti.  If  the  blades  and 
disk  are  not  cooled  then  Td,  and  are  the  same. 

These  simplified  algorithms  were  constructed  to  replace  the  original 
turbine  model  which  has  been  used  in  our  reference  system  models.  The 
original  model  does  an  excellent  job  of  estimating  the  stage-by-stage 
mass  and  performance  of  a  turbine,  and  it  was  a  significant  breakthrough 
for  our  modeling  effort.  But,  it  requires  multiple  iterations  to 
converge  on  the  proper  pressure  ratio  within  the  system  models,  and  it 
requires  substantial  computer  time.  The  simplified  models  avoid  all  of 
these  problems  while  using  the  original  model  as  a  basis.  Thus,  it  is 
fast  and  simple,  and  it  approximates  the  accuracy  of  the  original  basic 
phenomenological  model  over  limited  parameter  ranges . 

To  formulate  the  simplified  models,  I  assumed  that  the  logarithm  of 
turbine  mass  depends  on  the  logarithm  of  each  of  the  turbine  parameters 
independently,  and  that  the  dependence  on  each  is  linear,  except  for  the 
dependence  on  temperature.  That  is,  if  turbine  mass  is  plotted  as  a 
function  of  an  individual  parameter  on  a  log- log  plot,  the  result  is  a 
straight  line.  The  turbine's  mass  dependence  on  inlet  temperature  was 
assumed  to  be  exponential  when  the  blades  and  disk  are  not  cooled. 
Coefficients  for  each  of  the  parameters  were  found  by  making  the  best 
fit  to  the  mass  values  obtained  from  Hudson's  turbine  model. 


Using  this  procedure,  we  developed  simplified  mass  algorithms  for  the 
following  turbines: 

Nickel  --  H2-O2  combustion  product 

Carbon-Carbon  Composite  --  H2-O2  combustion  product 

Nickel  --  Hydrogen 

Carbon-Carbon  Composite  --  Hydrogen 

Nickel  --  Helium-Xenon 

Carbon-Carbon  Composite  --  Helium-Xenon 

Nickel  --  Potassium  vapor 
Tantalum  alloy  --  Potassium  Vapor 

These  aglorithms  were  fairly  good  over  limited  parameter  ranges,  and 
they  have  been  used  in  our  system  models.  However,  extreme  care  must  be 
taken  in  their  use.  Misleading  results  can  be  obtained  outside 
appropriate  parameter  ranges.  Even  within  appropriate  ranges,  problems 
can  arise:  blade  and  disk  strength  limitations  can  be  exceeded, 
unrealistic  blade  lengths  are  possible,  results  may  correspond  to  an 
excessive  number  of  stages,  and  there  are  other  potential  problems 
associated  with  fluid  velocity  and  blade  length  to  disk  diameter  ratios. 
We  always  confirm  system  results  with  the  original,  more  comprehensive 
model . 

Because  of  the  potential  problems  associated  with  using  the  simplified 
algorithms,  I  have  not  reproduced  them  in  this  report. 
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APPENDIX  D 


Cryogen  Storage 

The  cryogenic  storage  subsystem  consists  of  liquid  hydrogen  coolant  or 
liquid  oxygen  used  for  combustion,  a  tank,  multilayer  tank  insulation  a 
refrigeration  unit,  and  a  meteoroid  shield. 

Hydrogen  Weight 

The  mass  of  hydrogen  to  be  stored  is  equal  to  its  mass  flow  rate  through 
the  turbine  or  the  weapon,  whichever  is  greater,  multiplied  by  the 
system  operating  time.  The  mass  of  liquid  hydrogen  is  denoted  by  M. 


M  =  mT  .  (D-1) 

T  is  system  operating  time  in  seconds  and  in  is  the  hydrogen  flow  rate. 
The  volume  of  hydrogen  is  denoted  by  V. 


V  =  -  ■  .  (D-2) 

P 

p  is  the  density  of  the  stored  liquid.  The  hydrogen  is  stored  in  a 
spherical  aluminum  tank  at  20°  K  and  0.1  MPa  (1  Atm).  The  surface  area 
of  the  tank  is  A. 

A  =  47r  [  1^  ]  “  4.84V^/^  .  (D-3) 

Hydrogen  Tank  Weight 

The  tank  weight  can  be  found  using  a  stress  analysis.  The  stress  in  the 
tank  wall  cr(N/m2)  can  be  found  as  follows; 


wr2p 
^  2irrt 

where  r  is  tank  radius  (m) , 

P  is  pressure  (Pa) ,  and 
t  is  wall  thickness  (m) . 


(D-4) 


tank  weight 

tank  weight 
tank  volume 


47rr  tp 

1,5  Pp 
a 


(D-5) 

(D-6) 

(D-7) 


Values  for  these  parameters  for  an  aluminum  tank  follow.  The  tank  is 
aluminum  with  an  ultimate  strength  of  500  MPa  and  a  safety  factor 
of  2.5. 
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P  =  10^  Pa, 


p  =  2700  kg/m^, 
a  =  200  MPa,  and 


tank  weight 
tank  volume 


2.03  kg/m^. 


(D-8) 


Refrigeration  Unit  Weight 

The  refrigeration  unit  consists  of  refrigeration  equipment,  a  radiator, 
and  a  power  supply.  The  refrigeration  equipment  comprises  compressors, 
heat  exchangers,  turbines,  and  generators.  The  power  required  to  run 
the  refrigerator  system  is  found  by  dividing  the  heat  to  be  removed  by 
the  refrigerator's  COP.  We  will  assume  that  the  COP  is  15  percent  of 
the  Carnot  COP  for  cooling  loads  above  500  W  (Hudson,  1987). 

COP  =  0.15  =— TjT-  ,  (D-9) 

H  L 


P 


Q 

COP 


T  -T 
H  L 

0.15  T, 


(D-10) 


where  Tl  is  the  low  temperature  in  the  refrigerator  cycle  --  the 
temperature  at  which  the  cryogen  is  stored,  20  K, 

Th  is  the  high  temperature  in  the  cycle  --  the  radiator 
temperature,  in  K, 

Q  is  the  heat  that  must  be  removed,  in  W,  and 
P  is  the  power  required  by  the  refrigeration  unit  in  W/m^. 

The  total  heat  that  must  be  dissipated  is  equal  to  the  sum  of  P  and  Q. 
The  radiator's  area  can  be  found  by  dividing  the  heat  to  be  dissipated 
by  the  radiator's  emissive  power.  The  radiator's  weight  can  then  be 
found  by  multiplying  its  area  by  5  kg/m^,  the  radiator's  specific 
weight,  and  by  1.25,  a  factor  that  accounts  for  meteoroid  losses. 


radiator  6.25(P+Q) 

“  .88a(T^  -  250"^  ' 
I  ”  J 


(D-11) 


Here,  we  asstime  that  the  radiating  temperature  of  space  is  250  K,  and 
that  the  radiator's  emittance  is  0.88. 

We  calculate  the  refrigeration  unit's  power  source  weight  using  30  kg/kW 
which  is  the  specific  weight  of  a  proposed!  SP-100  power  system. 


D-2 


Thus,  the  power  source  weight  is  given  by  Equation  D-12. 

power  source  weight  -  0,03P  =  0.2Q(Th-Tl)/Tl  (D-12) 

<«- 

We  assume  that  the  refrigeration  equipment  weight  is  given  by 
equation  D-13  (Hudson,  1987),  This  weight  is  based  on  a  Garrett, 
reverse  Bray ton  refrigerator. 

refrigeration  -  -T  )/280  (D-13) 

equipment  weight 

The  total  refrigeration  unit's  weight  is  the  stun  of  the  radiator,  power  i 

source,  and  refrigeration  equipment's  weight. 

Weights  for  these  components  are  shown  in  Table  D-1  for  various  values 
of  Tg,  for  Tl  equal  to  20  K,  and  for  Q  values  of  100  and  1000  W. 

The  "optimum"  radiator  temperature,  Th,  is  345  K  at  100  W  and  355  K  at 
1000  W.  Thus,  the  optimum  temperature  is  not  very  dependent  on  the 
cooling  load,  Q,  and  I  selected  355  K  to  use  for  Th-  With  this  value  of 
Tg,  the  refrigeration  unit's  weight  is  given  by  equation  D-14. 

refrigeration  -  4.52Q  +  109xlOO  O*®®(J-°8(Q))**2.9  (D-14) 

unit  weight 

Insulation  Weight 

CRC  (1976)  gives  multilayer  insulation  conductivities  ranging  from  0.04 
to  0.2  mW/m-K  and  an  insulation  density  of  80  kg/m^.  The  heat  gain 
through  this  insulation  depends  on  the  difference  between  the  cryogen's 
temperature  (20  K)  and  the  temperature  of  space  (assumed  to  be  250  K) . 

(250-20)kA 

Q - ^  (D-15) 

t 

where  k  is  the  insulation's  conductivity  (we  used  0.0001  W/mK) , 

A  is  the  tank's  surface  area,  and 
t  is  the  insulation's  thickness  in  m. 

Insulation  weight  =  80  At  ,  (D-16) 

Insulation  thickness  affects  the  weight  of  the  refrigeration  unit  since 
it  determines  the  value  of  Q.  Thus,  to  optimize  insulation  thickness, 
we  must  minimize  the  combined  insulation  and  refrigeration  unit  weight. 

Table  D-2  shows  insulation  and  refrigeration  unit  weights  as  functions 
of  insulation  thickness  for  two  tank  sizes,  100  m^  and  1000  m^. 
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Table  D-1 

Refrigerator  Unit  Weights 


Radiator 

Temperature 

Tt,  (K^ 

Power 

Source 

mass 

Refrigeration 

Equipment 

mass 

Cooling  Load, 

Q  -  100  W 

300 

280 

204 

340 

320 

233 

345 

325 

236 

350 

330 

240 

355 

335 

244 

375 

355 

258 

400 

380 

276 

Cooling  Load, 

Q  -  1000  W 

300 

2800 

1235 

345 

3250 

1433 

350 

3300 

1455 

355 

3350 

1476 

360 

3400 

1500 

365 

3450 

1522 

375 

3550 

1566 

400 

3800 

1676 

Radiator 
mass  (ka) 


280 

142 

133 

125 

117 

94 

73 


2802 

1328 

1246 

1172 

1105 

1044 

937 

733 
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Total 
mass  (kg) 


764 

695 

694 

695 

696 
707 
729 


6837 

6011 

6001 

5998 

6005 

6016 

6053 

6209 


Table 
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Insulation  Thickness  Optimization 

Insulation 

Insulation 

Refrigeration 

Total 

thickness 

weight 

weight 

weight 

(ke) 

(ke) 

Tank  volume, 

V  =  100  m3 

.01 

83 

1492 

1575 

.02 

166 

808 

974 

.03 

350 

578 

828 

.035 

291 

512 

803 

.04 

333 

463 

796 

.045 

374 

424 

798 

.05 

416 

393 

809 

Tank  volxune, 

V  -  1000  m3 

.01 

386 

6679 

7065 

.02 

773 

3326 

4099 

.03 

1159 

2247 

3406 

.035 

1352 

1942 

3294 

.04 

1546 

1714 

3260 

.045 

1739 

1537 

3276 

An  insulation  thickness  of  0.04  m  minimizes  the  total  weight  for  both 
cases.  Using  this  insulation  thickness,  we  can  specify  the  cooling 
load  Q,  power  P,  refrigeration  unit  weight,  and  insulation  weight  in 
terms  of  tank  area  A. 


Q  *  0.575  A  in  watts 

• 

(D-17) 

P  «  64.2  A  in  watts 

(D-18) 

Insulation  weight  ='3.2 

A 

(D-19) 

Refrigeration  2.60A  + 
unit  weight 

109x10° .  0^68  (log  ( .  575A)  )**2 . 9 

(D-20) 

Meteoroid  Shield  Weight 

The  thickness  of  an  aluminum  meteoroid  shield  that  has  only  a  1  percent 
chance  of  being  penetrated  by  a  meteoroid  during  seven  years  in  low 
earth  orbit  is  specified  by  equation  D-21. 
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t  ■=  5.2  X  10-6(AT/-lnP)'»-29 


(D-21) 


where  T  is  the  exposure  time  in  seconds  (7  years), 

P  is  the  probability  of  no  meteoroid  penetrations  (.99), 
t  is  the  shield's  thickness  in  m,  and 
A  is  the  tank's  area  in  m^. 

This  comes  from  Fras  (1986)  using  parameters  for  an  aluminum  shield  and 
low  earth  orbit.  While  the  shield  will  stop  natural  meteoroids,  it  will 
not  protect  against  space  debris.  Aluminum  density  (2700  kg/m^)  is 
multiplied  by  tank  area  and  shield  thickness  to  get  shield  weight. 

shield  weight  -  14  A^-^^  (D-22) 


Hydrogen  Summary 

Using  equation  D-3,  the  weights  of  the  cryogen  storage  system  can  be 


summarized  as  follows: 

cooling  load  Q  -  2.78  V2/3  .  (D-23) 

refrigeration  power  =  311  V2/3  .  (D-24) 

tank  weight  =  2.03  V  .  (0-25) 

insulation  weight  -  15.5  V2/3  .  (D-26) 

refrigeration  unit  weight  = 

12.6  V2/3+109x100  0'*®8(^°8^2.78V**2/3))**2.9  .  (D-27) 

meteoroid  shield  weight  =  107  .  (D-28) 


V  is  the  tank's  volume  in  m®  and  all  weights  are  in  kg. 

The  oxygen  storage  subsystem's  mass  was  estimated  in  a  similar  manner  to 
that  used  for  hydrogen.  We  assumed  that  oxygen  is  stored  at  90  K  and 
0.1  MPa  in  a  spherical  aluminum  tank  surrounded  by  multifoil  insulation 
and  an  aluminum  meteoroid  shield.  We  also  assumed  that  the  oxygen  s 
refrigeration  system  will  be  integrated  with  the  hydrogen's 
refrigeration  system. 

Oxygen  Mass 

The  density  of  liquid  oxygen  at  0.1  MPa  is  1142  kg/m®. 

Oxygen  Tank  Mass 

The  algorithm  for  this  is  exactly  the  same  as  for  a  hydrogen  tank. 
Refrigeration  System  Mass 

The  oxygen  will  most  likely  not  have  a  separate  refrigeration 
system.  It  will  "piggy-back"  on  the  hydrogen's  refrigeration 
system.  We  will  assxune  that  the  added  refrigerator  mass  can  be 
estimated  as  follows: 
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X  H2  refrigerator  mass.  (D-29) 


Oo  refrigerator  mass  =  ^ 

^  ^  Qh  COPo 

Qo  is  the  heat  load  to  the  oxygen  tank. 

QO  *  kA(Ts  -  To)/t  ,  (D.30) 

where  k  is  insulation  conductivity  -  .0001  w/mK, 

A  is  oxygen  tank  area  - 

Vq  is  oxygen  tank  volume  in  m^, 

Tg  is  the  background  temperature  -  250K, 

Tq  is  the  oxygen  temperature  -  90K,  and 
t  is  insulation  thickness  in  m.  Qo  0 ,077V^  /t 
Qh  is  the  hydrogen  heat  load.  Qh  ==  2.8Vh®*^^^ 

Vh  is  the  hydrogen  tank's  volume. 

COPh  is  the  Carnot  COP  for  the  hydrogen  refrigerator. 

20 

COPh - .  (D-31) 

Th  -  20 

Th  is  the  refrigerator's  radiator  temperature,  355K. 

COPq  is  the  Carnot  COP  for  a  refrigerator  operating  between  Th  and  the 
oxygen  storage  temperature,  90K. 


90 

COPo - 

Th  -  90 

Thus,  refrigerator  mass  is  estimated  as  follows: 


(D-32) 


.0048 

Refrigerator  mass  =  -  x 

t 


-0.667 

Vo 

—  X  H2  refrigerator  mass  (D-33) 

Vh 


Insulation  Mass 

Multifoil  insulation  has  a  density  of  80  kg/m^. 

Insulation  mass  =  80  x  4.84Vo®-®®^  x  t  .  (D-34) 

To  find  the  best  insulation  thickness,  I  assumed  that  the  masses  of 
oxygen  and  hydrogen  are  equal.  (This  is  true  if  the  combustion  power 
generation  system's  turbine  inlet  temperature  is  around  1200K  which 
requires  the  combustion  of  roughly  equal  portions  of  H2  and  O2.) 

The  following  table  shows  how  refrigeration  system  mass  and  insulation 
mass  trade-off  for  two  oxygen  tank  sizes,  1.0  and  30.0  m^. 
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Meteoroid  Shield 

The  algorithm  for  meteoroid  shield  mass  is  exactly  the  same  as  for 
the  hydrogen  tank. 
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APPENDIX  E 


Hydrogen -Oxygen  Combustion 

The  heat  of  combustion  for  hydrogen  and  oxygen  is  13.4x103  KJ/kg  of 
(gas  reactants  and  gas  products) .  Hydrogen  is  assumed  to  enter  the 
combustion  chamber  at  temperature  and  it  requires  heat  to  be  heated 
from  Ti  to  the  reaction  temperature  which  is  300  K.  It  is  preheated 
when  it  cools  the  weapon,  power  conditioning  unit,  and  generator.  The 
oxygen,  on  the  other  hand,  must  use  some  of  the  combustion  energy  to 
heat  it  from  a  cryogenic  liquid  to  T^.  The  remaining  combustion  energy 
heats  the  combustion  products,  steam  and  hydrogen  (since  we  are  using 
excess  hydrogen  to  reduce  the  combustion  product  temperature)  from 
temperature  Tj  to  temperature  T.  The  energy  balance  equation  for  this 
process  is  given  by  equation  E-1. 


m'H[hH(T)-hH(Tr)]  +  ms[hs(T)-hs(Tr)] 

—  ms(13,400)  -  mo[ho(Tj.) -ho(Tso)  ]  -  mg[hjj(Tj.) -hH(Ti)  ]  +  moSp^  , 


(E-1) 


where  m'g  is  the  mass  of  noncombusted  hydrogen  in  kg, 
hg  is  the  entHalpy  of  hydrogen, 
mg  is  the  mass  of  steam  in  kg, 
hg  is  the  enthalpy  of  steam, 
mo  is  the  mass  of  oxygen  in  kg, 

ho  is  the  enthalpy  of  oxygen  being  heated  through  a 
supercritical  process, 

Spo  is  the  oxygen  pumping  energy  in  J/kg, 

Tgo  is  the  temperature  at  which  oxygen  is  stored  as  a  liquid 
-90  K,  and 

mg  is  the  mass  of  hydrogen  entering  the  combustion  chamber. 

We  want  to  find  the  mass  ratio.  Mg,  of  hydrogen  to  oxygen  that  results 
in  combustion  product  temperature  T,  the  turbine  inlet  temperature. 

In  the  combustion  process ,  0 . 125  kg  of  Hg  combines  with  1  kg  of  Og  to 
form  1.125  kg  of  steam.  Using  this,  we  can  rewrite  equation  E-1. 


(mH-0.125mo)[hg(T)-hg(Tr)]  +  1.125mo[hg(T)-hs(Tr)]  (E-2) 

=  1.125mo(13,400)-mg[hH(Tr)-hH(Ti)]-mo[ho(Tr)-ho(T3j]+moSp„  . 


mg  is  the  total  mass  of  hydrogen  entering  the  reaction.  We  can  divide 
both  sides  of  E-2  by  mg  to  get  Equation  E-3. 


„  =  0 ■  125  [hH(T)  -hH(T,)  ]+l .  125(13,400)  -ho(T,.)+hn(T„J  -1 . 125  [hg(T)  -h.(T.)  1+S.. 

hg(T)  -  hg(T,) 
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This  defines  the  needed  ratio  of  hydrogen- to- oxygen  depending  on 
hydrogen  inlet  temperature,  pump  power,  and  the  desired  combustion 
product  temperature . 

The  properties  of  the  combustion  products  depend  on  the  masses  of 
hydrogen  and  steam  in  the  combustion  products. 

m'H  ”  (^  •  0.125)mo 

ms  “  1.125  toq 

Combustion  product  mixture  properties  are  given  by  equations  E-6, 
and  8. 


h 


hsm's  +  hsmg  (MR-0.125)hH  +  1.125hs 

m'fl  +  nis  1  +  Mr 


(MR-0.125)gH  +  1.125gs 

7  -  - 

1  +  Mr 


(Mr-0.125)Rh  +  1.1-25Rs 

R  =  - - - - 

1  +  Mr 


7  is  the  specific  heat  ratio 

(1.4  for  hydrogen,  1.24  for  steam),  and 

R  is  the  gas  constant  in  kJ/kgK. 

(4130  J/kgK  for  hydrogen,  46  J/kgK  for  steam) 

Values  for  enthalpy  are  found  in  Appendix  G. 


(E-4) 

(E-5) 

f 

t 

(E-6) 

•(E-7) 

(E-8) 


E-2 


APPENDIX  F 


Radiator  Area 

A  radiator  consists  of  side-by-side  heat  pipes.  The  evaporator  end  (the 
end  where  heat  enters)  of  each  pipe  sticks  into  a  heat  exchanger  where 
it  absorbs  waste  heat  from  the  working  fluid  of  a  thermodynamic  power 
generation  cycle. 

Bravton  Radiator 

Working  fluid  from  the  turbine  transfers  heat  to  the  finned  evaporator 
ends  of  the  radiator's  heat  pipes.  Equations  F-1  and  F-2  describe  the 
energy  balance  in  the  heat  exchanger. 

dT 

-mC  -  hP  (T  -  T  )  ,  (F-1) 

P  dX  ^  ^  ® 

-  T^)  =  P£a(T^  -  T^)  .  (F-2) 

m  is  the  working  fluid  flow  rate, 

Cp  is  the  working  fluid's  specific  heat, 

Tf  is  the  working  fluid's  temperature, 

X  is  the  dimension  along  the  heat  exchanger, 
h  is  the  convection  coefficient  in  the  heat  exchanger, 

Px  is  the  heat  exchanger's  surface  area  per  unit  length  of  heat 
exchanger , 

Tg  is  the  surface  temperature  of  the  heat  pipe's  evaporator, 

P  is  the  radiator's  radiating  surface  area  per  unit  length  of  heat 
exchanger , 

£  is  the  radiator  emittance  (we  use  0.88), 
a  is  the  Stephan-Boltzmann  constant, 

T  is  the  radiator's  mean  sur'face  temperature,  and 
Tg  is  the  radiating  temperature  of  space  (we  use  250  K) . 

We  assume  that  the  temperature  drop  between  the  evaporator's  surface  and 
the  mean  radiating  temperature  is  20  K;  that  is,  Tg  -  T  =  20®K. 

With  this,  we  can  solve  equations  F-1  and  F-2  as  follows: 

hP3t(Tj  -  T  -  20)  -  P£a(T^  -  T^)  ,  (F-3) 


T^)  +  T  +  20  ,  (F-4) 

s 


F-1 


dTj  h?ca  2  dT  dT 
dX  hPx  dX  dX 


(F-5) 


r  4Pea  ]  dT 

■ihC  -  T  +  1 - P£a(T  -  T  ) 

P  [  hPx  J  dX 


(F-6) 


Letting  a  =  — ^ 


we  get  Equation  (F-7) 


Pdx  =  .  - -  +  - - ^ 

I  T  -  T  T  -  T 
s  s 


(F-7) 


The  integral  of  PdX  is  the  radiation  area,  A,  and  can  be  found  by 
integrating  the  right  side  of  Equation  (F-7). 


^  ,  M^out-^  V^^in  -  V  -  -1  "out 

A - 3  -  In  - .  T  )(f - ^  T 

4T  ^  out  s''^  in  s'^  J  [  s 

s 


1  I  T. 

-1  in 


T7  -  T 

in  •  s 


(F-8) 


Tout  mean  radiator  temperature  at  the  radiator's  fluid  outlet  end 

and  is  its  mean  temperature  at  the  inlet. 

These  values  can  be  found  by  iterating  to  solve  Equations  (F-9)  and 
(F-IO). 


^  A  A. 

T^.  =  T.  +  —  (TT  -  T  )  +  20 

fin  in  ,  „  in  s 

hPx 


(F-9) 


Zl  a. 

T;  =  t  +  -  (T  \  -  T  )  +  20 

fout  out  „  out  s 
hPx 


(F-10) 


Tfin  Tfout  inlet  and  outlet  working  fluid  temperature. 


F-2 


hPx 


Notice  that  the  term  P  is  always  grouped  together.  Px/P  is  the  ratio 

of  heat  exchanger  area  to  radiator  area.  Px/P  is  the  most  important 
factor  in  determining  the  difference  between  the  fluid's  temperature  and 
the  radiator  evaporator's  surface  temperature. 

The  above  analysis  was  first  formulated  by  A1  Juhasz  at  NASA  LeRC. 


Thermionic  Radiator 


This  analysis  is  similar  to  the  previous  one  but  it  assumes  that  the 
temperature  difference  between  the  reactor's  cooling  fluid  and  the 
radiator's  evaporator  surface  is  20 ®K  and  that  the  temperature 
difference  between  the  evaporator's  surface  and  the  mean  radiator 
temperature  is  20 ®K  for  a  total  difference  of  40 ®K. 


dT 

-mCp  =  P£a(T^  -  2501 


(F-11) 


Since  Tf  ==  T  +  40,  dTf 


dT,  equation  F-11  can  be  written  as  follows: 


P£Cr 


dT 


(F-12) 


Integrating  PdX  gives  us  the  radiator's  radiating  area,  A. 


A 


3 

4£aT 

s 


f  'T  ^+T  1 

r  T.  +T  1 

out  s 

1  -n 

in  s 

T  ^-T 

“  in 

T.  -T 

out  s 

in  s  ^ 

+  2  tan 


T 

'  T. 

out 

-  2  tan  ^ 

in 

T 

T 

s  J 

L  s  J 

(F-13) 


Tin  Tout  the  mean  radiating  temperatures  at  the  fluid  inlet  and 
outlet  ends  respectively. 


Tin  -  Tfi„  -  40 


(F-14) 


^out  —  ^fout  ■ 


(F-15) 


^fin  Tfout  fluid  inlet  and  outlet  temperatures  respectively. 
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Ranklne  Radiator 


We  assume  that  the  radiator  for  a  Ranklne  system  extracts  heat  from  an 
Isothermal,  condensing  fluid  and  rejects  It  to  space  at  an  Isothermal 
radiating  temperature.  The  area  for  this  radiator  Is  given  by 
equation  F-16. 

A - 2 - ^  .  (F-16> 

ea(T  -  T^) 

^  s 

Q  is  the  thermal  power  (watts)  that  must  be  dissipated. 

We  assume  that  the  temperature  drop  between  the  condensing  fluid  and  the 
evaporator's  surface  is  10  K  and  that  the  drop  between  the  evaporator 
and  mean  radiating  temperature  is  20  K,  thus  T  «  Tf  -  30. 
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APPENDIX  G 


Enthalpies  for  Hydrogen,  Oxygen,  and  Steam 


Hydrogen 

Jones  and  Hawkins  (1960)  give  the  following  expression  for  hydrogen's 
specific  heat: 

Cp  -  12.1  +•  .00218T  +  31.2/Jt  ,  in  kJ/kg-K  .  (G-1) 

Their  values  were  given  in  English  units  for  300  K  and  up.  From  this, 
we  can  find  enthalpy  by  integrating: 

h  -  12. IT  +  .001091^  +  62.4Jt  +  C  .  (G-2) 

At  low  temperatures  we  used  the  following  data  from  Reynolds  (1979)  for 
parahydrogen: 


T  (k) 

P  (’MPa') 

h  rKJ/ke^ 

20 

.1 

50 

100 

10 

1293 

200 

10 

2984 

300 

10 

4561 

We  used  these  because  most  of  our  processes  start  at  20  K  and  0.1  MPa, 
are  pumped  to  high  pressure,  then  heated  in  various  components.  To  the 
300  K  enthalpy  we  added  229  kJ/kg  to  account  for  75%  of  the  hydrogen 
changing  from  para  to  ortho.  We  assumed  the  transition  came  between  200 
and  300  K  although  we  realize  the  transition  depends  on  many  factors. 

We  fit  these  data  and  obtained  the  following  relation: 

h  -  .007T2  +  14. 7T  -  247  ,  in  kJ/kg.  (G-3) 

We  use  this  enthalpy  expression  below  300  K  and  Equation  (G-2)  above 
300  K.  To  make  them  continuous,  the  constant  term  in  Equation  (G-3)  is 
set  equal  to  -15, 

Because  enthalpy  is  rather  dependent  on  pressure  at  low  temperatures, 
the  accuracy  of  Equation  (G-3)  is  not  good  for  some  pressure  values  at 
temperatures  below  100  K. 

Oxygen 

For  oxygen  we  fit  data  from  Reynolds  (1979) .  Again,  we  assumed  that 
oxygen  would  start  at  0.1  MPa  and  be  ptimped  to  high  pressure  before 
being  used. 

h  -  1.0625T  +  117.5  ,  T  >  200  K  .  (G-4) 


G-1 


h  -  2.455T  -  161.3  ,  T  <  200  K  .  (G-5) 

These  are  in  kJ/kg. 


Steam 

For  steam,  we  used  the  expression  for  Cp  from  Jones  (1960)  and 
integrated  to  find  enthalpy. 

h  -  4.61T  -  206jT  +967  InT  -  787  (G-6) 

This  is  in  kJ/kg. 
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